SF  298  MASTER  COPY 


KEEP  THIS  COPY  FOR  REPRODUCTION  PURPOSES 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  NO.  0704-0188 


Public  reporting  burderi  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructbns,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comment  regarding  this  burden  estimates  or  any  other  aspect  of  this  collection 
of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway. 
Suite  1204,  Arlington.  VA  22202-4302.  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0701-0188),  Washinoton.  DC  20503 


9/14/97 


Technical  Report 


Heat  Transfer  in  Microchannels 


Xianming  Wang 
Timothy  A.  Ameel 
Robert  O.  Warrington 


I  aa  A 1 1  Mil  I  mi*  I  Ai  •]  11  asffifi  9SI1 


Louisiana  Tech  University 
P.O.  Box  7923  T.S. 
Ruston,  LA  71272 


■aaeiJitettizii^rcwiAreiifiifinigircwiTci  reva  ^  tA  a  i  ai  ai  asKf i  acti 


U.S.  ARMY  RESEARCH  OFFICE 
P.O.  BOX  12211 

RESEARCH  TRIANGLE  PARK,  NC  27709-221 1 


AGENCY  REPORT  NUMBER 


The  views,  opinions  and/or  findings  contained  in  this  report  are  those  of  the  author(s)  and  should  not  be  construed  as  an  official 
Department  of  the  Army  position,  policy,  or  decision,  unless  so  designated  by  other  documentation. 


Approved  for  public  release;  distribution  unlimited. 


riisMisriTMar/AC 


The  objective  of  this  research  was  to  develop  an  analytical  solution  to  the  heat  transfer  problem  in  microchannels  with  slip- 
flow  and  an  isothermal  wall,  a  heat  transfer  problem  for  gases  at  low  pressures  or  in  extremely  small  geometries,  and  to  verify 
this  solution  experimentally.  In  this  investigation,  an  analytical  expression  for  the  velocity  distribution  with  slip-flow  was 
obtained  which  involved  the  Knudsen  number  Kn  in  an  infinite  series  form.  Kn  for  extremely  small  channels  may  become  large 
enough  to  significantly  affect  the  velocity  distribution  and  consequently  affect  the  heat  transfer  properties.  A  mathematical 
model  of  temperature  distribution  was  established  by  combining  the  energy  and  momentum  equations. 

A  new  technique  for  evaluation  of  eigenvalues  for  the  solution  of  the  heat  transfer  problem  in  microchannels  was  developed 
from  the  method  of  Frobenius.  The  computational  results  show  that  the  method  is  effective.  The  local  and  average  Nusselt 
numbers  were  found  for  0.005  <  Kn  <  0.30  with  aspect  ratio  a  =  1,  2/3,  1/2,  1/4,  and  1/8.  Experiments  for  helium  flow  in  a 
rectangular  microchannel  of  dimension  1 17  jtm  x  24  jim  x  63.5  jim  and  a  microtube  with  inside  diameter  of  52  ^.m  were 
conducted. 
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ABSTRACT 


The  objective  of  this  research  was  to  develop  an  analytical  soluUon  to  the  heat 
transfer  problem  in  microchannels  with  slip-flow  -  a  heat  transfer  problem  for  gases  at 
low  pressures  or  in  extremely  small  geometries,  and  to  verify  this  solution 
experimentally.  In  this  investigation,  an  analytical  expression  for  the  velocity 
dtstnbuaon  with  slip-flow  was  obtained  which  involved  the  Knudsen  (Kn)  number  in  an 
mfinite  series  form.  The  result  showed  that  the  velocity  always  increased  as  the  Knudsen 
number  was  increased.  The  Knudsen  number  for  exriemely  small  channels  may  become 
large  enough  to  affect  significanUy  the  velocity  distribution  and  consequently  afi-ect  the 
heat  transfer  properties.  A  mathematical  model  of  temperature  distribution  was 
estabUshed  by  combining  the  energy  and  momentum  equations.  A  series  solution  was 

obtained.  Also,  expressions  for  the  local  and  overall  Nusselt  numbers  were  derived  in 
terms  of  the  Knudsen  number  and  Graetz  number. 

A  new  rechnique  for  evaluation  of  eigenvalues  for  the  solution  of  the  heat  transfer 
problem  in  microchannels  was  developed.  TTris  method  was  based  on  the  construction 

ofamarrix.  The  compurational  results  showed  that  rhe  method  was  effective.  Tlrelocal 

values  and  average  Nusselt  number  were  found  for  Kn  from  0.005  to  0.3  with  aspect  ratio 

a  =  I,  M,  1/2,  1/4  and  1/8.  Experimenrs  for  helium  through  a  microchannel  with 

dtmensrons  ofll7pmx24|rmx63.5mmanda  microtube  with  inside  diameter  of  52 
Jim  were  conducted. 
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NOMENCLATURE 


Aw 

tube  surface  area  [m^] 

M 

a 

aspect  ratio  h/b 

ap^  coefficient  in  Eq.  (4.2) 

NUoo 

bp,q{ 

'm,n)  coefficient  in  Eq.  (4.6) 

Cn 

coefficient  in  Eq.  (3.5) 

coefficient  in  Eq.  (4.11) 

Nux 

Dh 

hydraulic  diaraeter[m] 

P 

^m,n 

VnAin,n  in  Eq.  (4.13) 

q 

F 

specular  reflection  coefficient 

(Ur-Ui)/(Uw-Ui) 

Q 

G 

G{x,y)  in  Eq.  (3.3) 

r 

Gv 

volume  flow  rate 

Re 

Gz 

Graetz  number,  (RePr(D/L)) 

Rg 

hx 

local  convective  heat  transfer 

Pr 

coefficient  [W/m“-K] 

T 

he 

average  convective  heat  transfer 

Tb 

coefficient  [W/m^-K] 

Tl 

H 

characteristic  dimension 

(AT)ln, 

k 

thermal  conductivity  [W/m-K]; 

ratio  of  Nuoo(/:n)/Nuoo(0) 

u 

Kn 

Knudsen  number,  (A/D) 

Ui 

L 

length  of  ducts[m] 

Mach  number;  molecular 
weight 

fully-developed  Nusselt  num¬ 
ber  coefficient,  (hcD/k) 

Nusselt  number,  (hcD/k) 
local  Nusselt  number,  (hcD/k) 
fluid  pressure  [Pa] 
heat  flux  per  unit  wall  area 
tW/m2] 

heat  transfer  rate  [W] 
wMo 

Reynolds  number,  (QuD/p) 
ideal  gas  constant  [J/kg-K] 
Prandtl  number,  (v/a) 

T(r,x),  temperature  [K] 
bulk  temperature  [K] 
temperature  at  g=L  [K] 
log— mean— temperature 
difference  (LMTD)  [K] 
velocity  in  |  or  x  direction  [m/s] 
average  streamwise  velocity  of 
the  incident  molecules  [ra/s] 


Uf  average  strearawise  velocity  of 

the  reflected  molecules  [m/s] 
Uw  average  streamwise  velocity  of 

the  surface  [m/s] 

V  velocity  in  T)  or  y  direction  [m/s] 

w  velocity  in  g  or  z  direction  [m/s] 

Wo  average  velocity  [m/s] 

X  dimensionless  width,  (^Dh) 

X*  dimensionless  length,  (x/L) 

y  dimensionless  width,  (t]/D//) 

z  dimensionless  distance,  (g/Dff) 

Greek  Symbols 

a  fluid  thermal  diffusivity,  (k/gc) 

[m^/s] 

O/  velocity  eigenvalues 

^  (l+a)/2 

P'  (l+l/a)/2 

Y  ratio  of  specific  heats 

6p,q  coefficient  in  Eq.  (4.6) 

eigenfunction  Eq.  (5.8) 

^  width 

i)  width 

distance  along  duct 


function  in  Eq.  (4.1) 

A 

coefficient  in  Eq.  (4.1) 

X 

eigenvalue 

Kn 

mean  free  path  [m] 

V- 

dynamic  viscosity[kg/m  s] 

V 

kinematic  viscosity  [m^/s] 

Q 

density  [kg/m^] 

Jt 

3.141592654 

e 

^""^wl/fTo-Tw)  dimensionless 

temperature 

% 

(Tb~T\v)/(To— T^)  dimension¬ 
less  bulk  temperature 

^B,L 

dimensionless  fluid  bulk 

temperature  at  g=L 

^LN 

dimensionless  LMTD 

©(x,y) 

function  expressed  as  Eq.  (4.2) 

Subscripts 

0 

atz  =  0 

B 

bulk 

m 

average 

s 

slip— flow 

w 

wall 
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CHAPTER  1 


INTRODUCTION 

The  objective  of  this  research  is  to  develop  an  analytical  solution  of  the  heat  transfer 
problem  in  microchannels  with  slip-flow  —  a  heat  transfer  problem  for  gases  at  low 
pressures  or  in  extremely  small  ducts,  and  to  verify  the  analytical  model  experimentally. 
To  do  this,  the  velocity  profile  with  slip-flow  must  be  found  first.  A  mathematical  model 
of  the  temperature  distribution  in  slip  flow  should  be  established  by  combining  the  energy 
and  momentum  equations.  And  finaUy,  an  effective  technique  for  evaluation  of  the 
eigenvalues  for  the  series  solution  should  be  developed  and  heat  transfer  experiments 
with  microchannel/microtube  should  be  conducted. 

1.1  Heat  Transfer  Problem  in  Ducts 

By  the  end  of  the  last  century,  the  problem  of  forced  convection  heat  transfer  in  a 
circular  tube  in  laminar  flow  gained  interest  because  of  its  fundamental  importance  in 
physical  problems  such  as  the  analysis  and  design  of  heat  exchangers. 

The  Graetz  problem  is  a  simplified  case  of  the  problem  of  forced  convection  heat 
transfer  in  a  circular  tube  in  laminar  flow.  With  the  assumptions  of  steady  and 
incompressible  flow,  constant  fluid  properties,  no  ’’swiri”  component  of  velocity,  fuUy 
developed  velocity  profile,  and  negligible  energy  dissipation  effects,  Graetz  (1883) 
originally  solved  this  problem  analytically.  The  solution  by  Graetz  involved  an  infinite 
number  of  eigenvalues,  and  in  his  paper,  only  the  first  two  eigenvalues  were  evaluated. 
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Since  the  accuracy  of  the  Graetz  solution  mainly  depends  on  the  number  of 
eigenvalues,  it  is  extremely  important  to  obtain  more  eigenvalues,  as  Tribus  and  Klein 
(1953)  pointed  out.  For  seventy  years,  the  research  on  this  problem  focused  mainly  on 
finding  more  eigenvalues.  Abramowitz  (1953)  employed  a  fairly  rapidly  converging 
series  solution  of  the  Graetz  equation  in  making  the  calculation  and  found  the  lowest  five 
values  with  much  more  accuracy.  Sellars  et  al.  (1956)  extended  the  problem  to  include 
a  more  effective  approximation  technique  for  evaluation  of  the  eigenvalues  of  the 

problem;  they  could  get  any  number  of  eigenvalues  as  needed.  This  work  solved  the 
Graetz  problem  completely. 

Dennis  et  al.  (1959)  studied  the  case  in  rectangular  ducts.  By  employing 
homogeneous  linear  algebraic  equations  and  Rayleigh  quotient,  they  devoloped  a 
technique  for  evaluation  of  the  eigenvalues  for  the  analytical  solution  of  the  problem  of 
forced  convection  heat  transfer  in  rectangular  ducts  with  different  aspect  ratios.  The  same 
results  were  obtained  by  other  researchers  numerically  (Shah  and  London,  1978). 

1.2  Heat  Transfer  Problem  in  Ducts 
in  Slip-Flow 

Applications  of  microstructures  such  as  micro  heat  exchangers  have  led  to  increased 
interest  in  convection  heat  transfer  in  micro-geometries.  Some  experimental  work  has 
been  done ,  such  as  the  experimental  investigations  in  microtubes  (Choi  et  al.,  1991),  in 
microchannels  (Pfahler  et  al.,  1991),  and  in  micro  heat  pipes  (Petersen  et  al.,  1993). 
Appropriate  models  are  needed  to  explain  the  significant  departures  in  the  micro-scale 
expenmental  results  from  the  thermofluid  correlations  used  for  conventional-sized 
geometries.  For  example,  Choi  et  al.  (1991)  conducted  heat  transfer  experiments  using 
essentially  smooth  tubes  with  a  relative  roughness  of 0.0003  and  a  diameter  ranging  from 
3  pm  to  8 1  pm.  As  shown  in  Fig.  1 . 1 ,  the  measured  heat  transfer  coefficients  in  laminar 
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flow  in  small  tubes  exhibited  a  Reynolds  number  dependence,  in  contrast  to  the 
conventional  prediction  for  fully  estabHshed  laminar  flow,  in  which  the  Nusselt  number 
is  constant.  Also,  an  experimental  investigation  of  fluid  flow  in  extremely  small  channels 
showed  that  there  are  deviations  between  the  Navier-Stokes  predictions  and  the 
experimental  observations  (Pfahler  et  al.,  1 99 1 ). 

Therefore,  some  effects  and  conditions  that  are  normally  neglected  when 
considering  macro-scale  flow  must  be  taken  into  consideration  in  micro-scale 
convection.  One  of  these  conditions  is  slip-flow  (Flik  et  al.,  1992,  Beskok  and 
Kamiadakis,  1992).  It  has  been  found  that  the  analytical  model  combined  with  slip-flow 
condiuons  can  fit  the  experimental  data  in  microchannels  with  a  uniform  cross-sectional 
area  (Arkilic  et  al.,  1994)  and  with  a  non-uniform  cross-sectional  area  (Liu  et  al.,  1 995). 

Shp-flow  occurs  when  gases  are  at  low  pressures  or  for  flow  in  extremely  small 
passages.  At  low  pressures,  with  correspondingly  low  densities,  the  molecular  mean  free 
path,  which  can  be  expressed  as  Eq.  (1.1),  becomes  comparable  with  the  body 
dimensions,  and  then  the  effect  of  molecular  structure  becomes  a  factor  in  flow  and  heat 
transfer  mechanisms  (Eckert  and  Drake,  1972). 

A  =  l^^uT 

"  Pyl2gcM  (1.1) 

The  relative  importance  of  effects  due  to  the  rarefaction  of  a  gas  can  be  indicated  by 
the  Knudsen  number,  a  ratio  of  the  magnitude  of  the  mean  free  molecular  path  in  the  gas 
to  the  characteristic  dimension  in  the  flow  field.  The  effects  of  rarefaction  phenomena  on 
flow  and  heat  transfer  becomes  important  when  the  Knudsen  number  can  no  longer  be 
neglected.  The  Knudsen  number  may  be  defined  as 

Kn  =  —  =  P  I^PuT 
P>H  P^ni^gcM 


(1.2) 
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In  defining  when  sUp-flow  occurs,  Beskok  and  Kamiadakis  (1992)  have  proposed 
to  classify  four  flow  regimes  for  gases,  as  follows: 


Continuum  flow: 
Slip-flow: 
Transition  flow: 
Free  molecular  flow 


Kn  <  10-3 

10-3  <^02<  0.1 

0.1  <Kn<\{} 

10  ^Kn 


When  slip-flow  occurs,  the  gas  adjacent  to  the  surface,  in  contrast  to  its  behavior  in 
continuum  flow,  no  longer  reaches  the  velocity  or  temperature  of  the  surface.  In 
continuum  flow,  intermolecular  colhsions  dominate  the  flow  field,  and  a  usual  boundaiy 
condition  (continuous  boundaiy)  at  the  interface  between  a  fluid  and  a  solid  surface  is  that 


the  fluid  adjacent  to  the  surface  assumes  both  the  velocity  and  temperature  of  the  surface. 

In  the  case  of  slip-flow,  the  molecular  mean  free  path  is  rather  larger  than  any 

significant  body  dimension  so  that  most  of  the  gas  molecules  striking  and  leaving  the 

body  surface  do  not  collide  with  free-stream  molecules  until  veiy  far  from  the  surface. 

Thus,  the  gas  at  the  surface  has  a  tangential  velocity,  and  it  appears  to  slip  along  the 
surface. 

The  slip  velocity  can  be  expressed  as  follows  as  a  function  of  the  velocity  gradient 
near  the  wall: 


“  %  V=o 


and  Arkilic  et  al.  (1994)  give  the  expression  as  follows: 


(1.3) 


or 


du/c 

dy/H 


(1.4a) 


u  X 

S  p 


(  —  ) 


(I.4b) 
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which  includes  the  consideration  of  three  accommodation  coefficients  represented  by  the 
speculation  reflection  coefficient  F,  which  has  values  that  typically  lie  between  0.9  and 
1  (Ebert  etal.,  1965).  In  the  case  of  F  having  the  value  one,  Eq.  (1.4)  becomes  Eq.  (1.3). 
For  simplicity  in  this  investigation,  Eq.  (1.3)  was  applied  to  evaluate  the  velocity. 

The  temperature  boundary  condition  can  be  regarded  as  discontinuous;  that  is,  there 
is  a  jump  in  temperature  at  the  wall.  Actually,  the  temperature  of  the  gas  near  the  solid 
surface  changes  continuously  to  the  temperature  of  the  surface  but  only  in  a  very  thin  layer 
on  a  microscale  so  that  on  the  macroscale  there  appears  to  be  a  jump  in  temperature 

between  the  surface  and  the  adjacent  gas.  Eckert  and  Drake  (1972)  give  expressions  for 
the  temperature  jump  condition: 


where  represents  the  mean  free  path  for  collisions  between  a  moving  molecule  and 
the  fixed  molecules  (or  is  the  average  number  of  colUsions  per  unit  distance  ) 
(Present,  1 958),  and  aT/dy  is  the  temperature  gradient  at  the  waU.  As  shown  in  Eq.  (1.5), 
the  temperature  jump  is  proportional  to  A^^,:  a  small  value  of  A^^  means  that  a  great 
numbers  of  molecules  are  involved  in  energy  transport  so  that  the  temperature  jump  is 
small,  while  a  large  value  of  A„  means  that  fewer  molecules  are  involved  in  energy 
transport  so  that  the  temperature  jump  is  relatively  larger.  When  A^  is  smaU  enough,  as 
in  a  conventional  case,  the  temperature  jump  may  be  neglected. 

With  the  introduction  of  the  slip-flow  condition  into  the  Graetz  problem,  it  becomes 
more  difficuit  to  solve  such  problems.  The  classical  Gmetz  problem  is  governed  by  a 
partial  differential  equation  with  a  continuous  temperature  boundaiy  condition;  while  in 
the  Graetz  problem  combined  with  slip-flow  condition,  the  temperature  boundaiy 
condition  is  no  longer  continuous,  which  makes  the  solution  as  well  as  the  cotresponding 
eigenvalues  much  more  complicated  and  difficult. 
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In  the  case  of  circular  tubes  in  slip-flow,  that  is,  gases  at  low  pressures  or  in  extremely 
small  tubes,  the  heat  transfer  coefficient  depends  not  only  on  the  Reynolds  number  and 
Prandtl  number,  but  also  on  the  Knudsen  number.  Barron  et  al.  (1995)  developed  a 
technique  and  evaluated  the  eigenvalues  of  the  analytical  solution  for  this  problem  in  the 
case  of  a  constant  wall  temperature.  Wang  et  al.  (1995)  solved  this  problem  completely. 
Their  studies  shows  that  the  Nusselt  number  increases  significantly  with  the  increase  of 
the  Knudsen  number,  as  shown  in  Fig.  1.2.  Ameel  et  al.  (1996)  studied  the  case  with 


Fig.  1.2  Fully  developed  Nuasa  function  of  [Barron  et  al.] 

(constant  wall  temperature) 


constant  heat  flux,  and  the  results  indicate  that  the  fuUy-developed  Nusselt  number 
decreases  with  the  increase  of  the  Knudsen  number,  as  shown  in  Fig.  1.3.  Tlierefore, 

some  research  is  needed  to  analyze  this  type  of  problem  in  the  case  of  rectangular  ducts 
in  slip>-flow. 
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1.3  Related  Research 

Ebert  and  Sparrow  (1965)  have  found  the  fluid  velocity  distribution  for  slip-flow  in 
microchannels.  Their  results  can  be  summarized  as  follows: 

The  momentum  equation  can  be  written  as  foUows  in  the  coordinate  system  shown 
in  Fig.  1.4: 

d^w  ,  d^w  _  1  dp 

dr]^  Mdg  (1.6) 


Fig.  1.4  Coordinate  system 
with  the  slip— flow  boundary  condition 

=  atr/  =  h,  Os|<i>  (1.7a) 


-^  =  0  at  rj  =  0,  0  ^  ^  <  b 

rv  f-  r.  ^ 

0  ar  ^  —  0,  0  ^  ij  <  h 


(1.7b) 

(1.7c) 


(1.7d) 


The  velocity  can  expressed  as 
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-  =  n  y  <^os  g t]  ^  sin  a: 

i=i  a?  ^\  +  2  Kn  sin^ 
f*  dz  ' 


cosh 


cosh  -^  +  2  Kn  Qi  sin  ^ 


and  the  mean  velocity  can  be  determined  by 


— .0  =  _0  y  fl  , _  O/  ,  Q; 

«• = 1  1  +2Kn  sin^ 

A*  dz 


sin^  O; 


cosh  %  +  2  Kn  sin  ~ 


where  the  eigenvalues  Oi  satisfy 


a  tan  a,-  =  — 

'  '  2  Kn 


(1.10) 


Dennis  e.  al.  (1959)  solved  .he  problem  of  forced  heat  convecdon  in  laminar  flow 
.hrongh  rectangular  ducts  wid.  non-slip-flow.  Their  msuiu  may  he  summarized  as 


follows: 


The  energy  equation  can  be  written  as  follows  in  the  coordinate 


system  shown  in  Fi^. 


/t(  _  QT 


(1.11) 


Fig.  1.5  New  coordinate  system  for  dimensionless  variables 


with  the  boundary  condition 
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T-h,  =  Tj  g>  0 

T  =  To  g<  0 


Introducing  dimensionless  variables 


X 


T>h 


y 


-  =  J_ 

^  Dfj 


d  = 


T-T^ 

To-T, 


(1.12a) 

(1.12b) 


where  the  hydrauHc  diameter  D„  =  4A/S,  and  the  aspect  ratio  n  and  Knudsen  number 


a 


Kn  = 


2-F^m 
F  2h 


When  the  aforementioned  change  of  variables  is  carried  out  and  neglecting  the  axial 
conduction  kd'^TId<;^  for  large  enough  Peclet  number,  since  it  is  of  order  (1/Pe)2 

compared  to  the  axial  convection  term  H-dT/Ss,  the  governing  equations  and  the  boundary 
conditions  become 


d^w  _  T>ffdp 
dx^  dy^  ~  M  dz 


or 


Vj^w 


F  dz 


wDh^ 
^  dz 


or 


_  w  BO 
^odz 


where  _  d2/Q^2  ^  ^2/^2  velocity,  and 


w  = 

-KnjS^ 

at 

y  =  0  and  y?. 

o 

A 

A 

w  = 

-Kn  y5'-^ 

at 

X  =  0  and  y9'. 

V 

V 

o 

dw  ^ 

=  0  at  y  = 

P 

dy  " 

2  ’ 

0  <  X  <  yS' 

dw  _ 

II 

O 

II 

dx 

2  ’ 

VI 

VI 

o 

(1.13) 

(1.14) 

(1.15a) 

(1.15b) 

(1.15c) 

(1.15d) 


12 


where 
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6g  - 


w  6  dx  dy 


r  r 

J  J  D' 


w  dx  dy 


(1.20) 


=  ZCn^e-^’'^ 

n=\ 


and  Nuoo  =  Ai  /4. 


Eigenvalues  An  can  be  determined  by  the  following  infinite  set  of  homogeneous 
linear  algebraic  equations 


w  oo 

J?o  ^  ^  (  m,/j  =  0,1,2,  ...  )  (1.21) 


p=0  q 

where 


dp,g{m,n)  =11  0m,n0p.9  dx  dy 


(1.22) 


^m,n(^>n)A;n,n  —  ( 1 4a )  ( nfi  +  a^rfi). 


(1.23) 


bp,q{m,n)  |„_^l  ^f"+p.n  +  q  ~  ^m  +  p,\n-q\  } 


(1.24) 


ffi  fP' 

w  —  w  .ijr;cx 

u,-,-  —  I  I  — -  rnx  r/ic  r — £. 


(-0-)  COS  i-jr)  dx  dy  (ij=0,1.2, ... )  (1.25) 


bp,q(m,n)  { 


^  0  for  all  p,q,m,n  0 
=  0  for  any  p,q,m,n  =  0 


(1.26) 
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Dennis  et  al.  (1959)  estimated  the  eigenvalues  for  different  aspect  ratios  for 
nonslip-flow  by  employing  the  Rayleigh  quotient 


fP'  , 

-  ^nVj  dndx  dy/ 

.  0  J  0 


rdn^dx  dy 


/>  =  0 


q=0 


p,q(^) 

(1.27) 


for  example,  for  aspect  ratio  a  =  I,  A]  =11.91  and  Nu«  =  2.98. 

Therefore,  based  on  those  research,  a  model  for  the  case  of  rectangular  channel  in 
slip-flow  can  be  established  and  a  technique  of  evaluation  of  the  corresponding 
eigenvalues  for  the  analytical  solution  of  the  model  should  be  developed. 


CHAPTER  2 


VELOCITY  AND  TEMPERATURE  DISTRIBUTIONS 


In  order  to  buUd  the  mathematical  model  for  the  problem  in  slip-flow,  the  velocity 
profile  must  be  found  first  in  the  new  coordinate  system.  In  this  chapter,  based  on  some 
assumptions,  the  expression  for  velocity  will  be  derived  from  the  continuity  equation  and 
momentum  equation.  The  slip  condition  wiU  be  used  to  evaluate  the  slip  velocity,  and 
the  velocity  will  be  expressed  in  tenns  of  a  Knudsen  number.  A  mathematical  model  of 
temperature  distribution  in  slii>-flow  will  be  established  by  combining  the  energy  and 
momentum  equations. 

2.1  Velocity  Distribution 

Consider  the  flow  of  a  fluid  in  a  rectangular  duct,  as  shown  in  Figure  2.1: 


L 

i 

1 

—  5  - 

- 1 - 

1 

— ^ — 
• 

T 

1. 

- — ib--H 

Fig.  2.1  Coordinate  system  for  the  problem 

For  this  model,  the  following  conditions  have  been  assumed  (Barron,  1996): 

(1)  The  flow  is  steady.  This  means  that  the  properties  of  the  flow  are  time 
independent. 

(2)  The  fluid  is  incompressible  (or,  if  a  gas  is  considered,  the  Mach  number  is  less 
than  0.30).  In  this  case,  the  density  may  be  assumed  constant. 
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(3)  The  flow  is  fully  established.  In  this  case,  the  axial  velocity,  w,  is  a  function  of 

the  coordinates  7  and  ^  only,  and  not  a  function  of  the  axial  coordinate  g.  In  addition,  the 
radial  velocity  is  zero. 

(4)  The  ’’swirl”  component  of  velocity  is  identically  zero.  This  means  that  both  u  = 
0  and  V  =  0. 

(5)  Fluid  properties  are  constant. 

(6)  Energy  dissipation  effects  are  negligible. 

(7)  The  tube  wall  temperature  is  constant 

2.1.1  Continuity  Equation 

The  general  continuity  equation  can  be  written  in  the  cartesian  coordinates  as 
follows: 


^  ^  =  0  (2.1) 

For  steady  flow  of  an  incompressible  fluid,  the  equation  above  reduces  to: 

—  +  —  4.  —  n 

^  (2.2) 

For  fiilly  developed  flow. 


Therefore, 


du  .  dv  _  f. 


Since  the  nomal  velocity  is  zero  at  the  walls  (the  wall  is  impermeable),  we  must  conclude 


that: 
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and 


w  =  0  (identically). 


V  =  0  (identically). 


2.1.2  Velocity  Distribution  with 
Slip  Condition 

The  Momentum  Equation  can  be  written  in  the  coordinates  shown  in  Fig.  2.2  as 
follows  (Kays  et  al.  1993)  for  fully-developed  (hydrodynamic)  flow: 

d^w  ^  d'^w  _  1  dp 

Qjjl  n  ^2.3) 


Fig.  2.2  New  coordinate  system  for  dimensionless  variables 


with  the  slip— flow  boundary  condition 

^  at  >7  =  0,2/1,  0  <  ^  < 

^  ^  -  0>2b,  0  <  T}  <  2h 

■^-0  at  T}  =  h,  0  ^  ^  <  2b 

dw  A  f. 

-  0  at  ^  =  b,  0  ^  T}  ^  2h 


(2.4a) 

(2.4b) 

(2.4c) 


(2.4d) 


w  = 

-Kn^lE 

dy 

at 

Vi 

II 

o 

o. 

o 

A 

V 

A 

(2.6a) 

w  = 

-KnP’^ 

dx 

at 

X  =  0  and  P',  0  <  y  < 

(2.6b) 

dw  _ 

=  0  at  y  = 

VI 

H 

VI 

o 

dy  ■ 

2  ’ 

(2.6c) 

dw  _ 

=  0  at  X  = 

VI 

VI 

o 

dx 

2  ’ 

(2.6d) 

The  solution  may  be  proposed  as 


w 

2jtL^ 

dz 


00 

=  E  w,  (x)  cos  (~-\)a-, 

/=  1  P 


(2.7) 


where  Oi  are  a  set  of  eigenvalues,  and  the  W,  (tc)  are  a  set  of  j  -  dependent  functions.  Eq. 
(2.7)  identically  satisfies  boundary  condition  (2.6c),  that  is 
dw 

Sy  _  2  V  2y 

55  ~  ^  aty  =  PI2 

dz 


By  substituting  Eq.  (2.7)  into  the  boundary  condition  (2.6a),  we  have 
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yW  dz 

-Kn  ^  ^  =  -Kn  (x)  sin  (^-l)a,  ] 

SO 

00 

IWi  (x)  [cos  i^-Da,  -  IKn  a,-  sin  (-|^-l)a,.  ]  =  0 
For  a  non-trivial  velocity  solution,  it  is  necessaiy  that 


a,  ,an  a,  =  ^ 

This  is  the  condition  from  which  the  eigenvalues  are  determined. 

Now,  we  find  the  function  W;  (x).  To  do  this,  Eq.  (2.5)  can  be  rewritten  a< 


d^w  I  d'^w  _ 
Bx^  By^ 


«  -  (f -»«,■ 


From  Eq.  (2.7),  we  have 


=  _£r±  y 
dx^  dzii,  dx^ 


2y 

cos  (-^-1)0; 


B^w 

By^ 


Diidp^-^ 

P  dz 


Sw,  cos  (%-l)a,. 

<=1  P 


SO  we  derive  a  differential  equation 


(fiW, 

dx^ 


2a  i  - 

-  (y)2  Wr  Qi  =  0 


(2.8) 


(2.9) 


(2.10) 


(2.11) 


(2.12) 
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We  can  obtain  the  homogeneous  solution  from 


=0 


that  is. 


(2.13) 


1  .2a,- 

A  =±(y) 

Therefore,  let 


(2.14) 


C,e-f^  +  C^e-1^^  + 

substituting  Eq.  (2.15)  into  Eq.  (2.12),  we  have 


(2.15) 


2a,-  - 

-  (■y)^C3  -Qi  =0 


C,  = 


(2.16) 


From  boundary  condition  (2.6d),  that  is. 


(2.17) 


Therefore,  we  obtain 


Cj  —  Cj< 


C,e7  [  ]  +  C, 


—  Dj  •^(2x— /S')  4-  C3 


(2.18) 
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And  from  boundary  condition  (2.6b),  we  have 


or 


w 


M  dz 


dw 

dx 


=  (x)  cos  (-^-1)0,-  =  -Kn  /5 

_  a>rV  •  .2y  ,,  , 


dW  (x) 

W;  W  +  Kn  =  0 


Substituting  Eq.  (2.18)  into  Eq.  (2.19),  we  have 


Cj[e ]  +  C^  +  c^Kn 


y5 


=  C^/[i€‘  +  \  )  +  2  Kn  ^(e^-1)] 


•>-£<>  «  ]  =  0 


4a, 2  ^ 


[  cosh^  +  2  Kn  ^sinh^  ] 


a,  .  ,  a, 

*  tf^v  mXa _ * 


Evaluation  of  i^-fx)  can  be  made  from  the  assumption 


(x)  cos  (-—-1)0:  =  1 

j=i  p 


and  the  property  of  orthogonality 


f  [cos  Oi  (^-1)]^  =  f 

JP/2  P  P  J  ^ 


^  . 2y  2y 

cos  a  I  (-^-1)  J(-^-l) 
^/2  P  P 


(2.19) 


(2.20) 


(2.21) 


(2.22) 
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Let  o  =  (2y/p  -  1),  we  have 


Zq  i  (x)  cos  a,  cr  =  1 


(2.23) 


From  the  orthogonality  we  have 


ij  I  cos  a  I  a  cos  a,-  o  da 

J  0 


cos  Oj  a  da 


(2.24) 


Carrying  out  the  integration  for  a[  =  aj ,  we  obtain: 
^  ^  a,-  +  sina,-  cosa,-  ,  sina. 

^  J - 7T— 


Q .  =  _  _2  ,  sinui- 

a,-  +  sinOi  cosa;  ^  2  Kn  sin^a,  ^ 


(2.25) 


Now,  the  velocity  can  be  expressed  as 


uJr  =  y  Qj  WM)  ^  sin  ai  cosh  ^(2x/^'-l) 

i  +  2  Kn  sirfi  a.-  Tn^h~^^~Z~n~r  ^ 

«  cosh  -^+  2  Kn  a,-  sin 

(2.26) 

where  C  =  (DH//x)(dp/dz)  and  the  mean  velocity  is  found  by  integration  of  Eq.  (2.26) 
across  the  section  of  the  channel:  thus, 


■il 


*^o  =  T  vv  flL4 


(2.27) 


Upon  carrying  out  the  indicated  operation,  we  obtain: 
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j/c  =  -2S  -2.(. 


sin^  a, 


O-i 

-)(7r- 


sinh.  ^ 


i-ia,-  co,h^^2Kna,,m^ 


(2.28) 


2.2  Temperature  Distribution 

The  general  temperature  field  equation  for  flow  of  an  incompressible  fluid  with  zero 
swirl  or  angular  components,  zero  energy  generation,  and  negUgible  frictional  energy 
dissipation  is: 


with  the  boundary  condition 


(2.30) 


Tu,  =  r, 


T  =  Tr. 


g  >  0 


?  <  0 


(2.31a) 

(2.31b) 


introducing  dimensionless  variables 


^  =  tI-  )■  =  -^  2  =  fl  =  I£L 

Du  D,,  T _ T 


To-T, 


When  the  aforementioned  change  of  variables  is  carried  out  and  the  axial  conduction 

can  bn  neglected  for  a  laige  enough  Peclet  number,  the  governing  equations  and 
the  boundary  conditions  become 


k(  ^4.  ^  dT 

av2  or 


V,^d  = 

^  ^Odz 


(2.32) 


where  V]2  =dW  +  dVd}’^,  with  the  boundary  condition 
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—  0  2  >  0 


(2.33a) 


0=1  z  <  0 


(2.33b) 


2.3  Summary 

In  this  chapter,  the  velocity  distribution  with  slip-flow  has  been  obtained  which  can 
be  expressed  simply  in  terms  of  the  Knudsen  number.  From  the  relationship  of  Kn  in 
terms  of  pressure  and  dimension,  we  can  see  that  Kn  in  microtubes  may  become  large 
enough  to  significantly  affect  the  velocity  distribution  and  consequently  affect  the  heat 
transfer  for  this  problem.  Also,  a  mathematical  model  of  temperature  distribution  in 
slip-flow  has  been  estabUshed  by  combining  the  energy  and  momentum  equations. 


CHAPTER  3 


ANALYTICAL  SOLUTION 


In  the  last  chapter,  the  velocity  distribution  was  expressed  in  terms  of  mean  velocity 
and  Knudsen  number,  and  a  mathematical  model  of  temperature  distribution  in  slip-flow 
was  established  by  combining  the  energy  and  momentum  equations.  In  this  chapter,  a 
series  solution  will  be  obtained  by  the  method  of  Frobenius.  Considering  the  boundai^ 
condition,  a  temperature  distribution  in  terms  of  a  generalized  Fourier  series  will  be 
denved.  Also,  expressions  for  the  local  and  overall  Nusselt  numbers  will  be  obtained. 

3.1  Graetz  Solution 

3.1.1  Separation  of  Variables  Solution 
The  governing  equation 

^odz  (3.1) 


e(x,y,z)  =  G(x,y)  Z(z) 
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Making  this  substitution  into  Eq.  (3.1)  and  rearranging  results  in  the  following: 


^  d;c2  dy2  Zdz~~'' 


where  A  is  an  arbitrary  constant.  The  ordinary  differential  equations  which  result  are: 


+  AZ=  0 


V.^G  +  ^ig^C=  0 


with  boundary  conditions: 


G(0,0)  =  0 


G(0,l)  =  0 


G(r,0)  =  0 


G(l’,l)  =  0. 


The  solution  of  Eq.  (3.4)  is: 


Z{z)  =  C  exp[  -Az  ] 


The  constant  C  in  Eq.  (3.6)  will  be  evaluated  below  and  A  will  be  evaluated  later. 


3.1.2  Determination  of  Constants  C„ 

We  can  wnte  the  solution  for  the  temperature  distribution  in  terms  of  a  generalized 
Fourier  series,  as  follows: 
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e  =  I  C„G„(x,y)e-^'^  .35) 

«=1  '  *  ^ 

where  G„  and  A„  are  eigenfunctions  and  eigenvalues  required  to  make  the  solution  to  the 

Eq.  (3.3)  subject  to  the  boundary  conditions  deduced  from  the  boundary  condition  Eq. 

(3.2a) 

G(O.O)  =  G(0,l)  =  G(r,0)  =  G(l’,l)  =  0. 

The  constants  Q  can  be  found  from  the  entrance  condition,  Eq.  (3.2b) 
atz  =  0,  q( x,y,0)=  1 


Making  this  substitution  into  Eq.  (3.5),  we  obtain: 

00 

^  (3.6) 

The  governing  differential  equation,  along  with  the  boundaiy  conditions,  is  a 
Sturra-Liouville  problem,  with  a  weight  function, 

w  =  w(x,y) 

where  w(x,y)  is  the  z-component  of  velocity  and  from  orthogonality. 


H'(A:,>0GnGm  dx  dy 


0 


for 


(3.7) 


The  constants  may  be  evaluated  by  multiplying  Eq.  (3.6)  on  both  sides  by 


^ix,y)G„^ 
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and  integrating  between  on  the  domain  D’  bounded  by  C’.  Only  the  term  in  which  m 

=  n  is  the  result  non— zero,  and  we  find; 


therefore. 


w{x,y){Gn)‘^dx  dy  = 


G„  dx  dy 


C„  = 


w{x,y)  G„  dx  dy 


M'U.y)  iG„)^dx  dy 


Each  C„  has  arbitrary  amplitude  which  we  choose  for  convenience. 


C„  = 


'  f 

w{x,  y) 

J  J  D' 

w{x, 

J  J  D’ 


y)  Gn  dx  dy 


(3.10) 


y)  dx  dy 


The  temperature  e(x,y,z)\s  therefore  known  to  any  desired  accuracy  once  sufficient 
c„  have  been  found.  Two  further  thennal  quantities  am  of  interest.  Experimental 
measurements  are  made  on  the  basis  of  a  mean  mixed  temperature  of  the  fluid,  that  is, 

efx,y,z)averagedwithrespect  to  the  local  fluid  velocity  over  any  section  of  thechannel. 

This  temperature  is  a  function  of  r  only  and  its  difference  between  any  two  sections  gives 
a  measure  of  the  heat  transferred  across  the  wall  between  them. 

From  these  expressions,  we  can  see  that  the  coefficients  C„  and  must  be 
predetermined  in  order  to  calculate  the  temperature. 


29 


30 


where  s  is  the  distajice  measured  along  the  perimeter  of  the  boundary  curve  C  in  an 

anti-elockwise  direction  and  v  is  the  outward-drawn  coordinate  normal  to  the  channel 

waU;  v-=v/Dh.  s  '=s/D„  .  Equating  the  heat  flux  from  Eqs.  (3.12)  and  (3. 1 3)  and  malting 

the  subsdmdons  from  Eq.  (3. 1 1)  for  the  temperature  gradient  at  the  wall  and  for  the  local 

bulk  temperature,  the  Mowing  expression  is  obtained  for  the  local  or  ••point”  Nusseit 
number. 

mz)  =  M  as'/se,  (3 ,,, 

We  now  eliminate  q  in  tenn  of  the  G.  by  Eq.  (3.5)  and  apply  Green's  theorem  to  Eq.  (3.3), 
SO  that 


^  ds'  = 
dv'  ^  Wo 


w(j:,y)  G„  dx  dy 


(3.15) 


U^ins  Eq.  (3.10)  and  since  //  p-  w  dx  dy  =  d  •»«,  where  d'  is  the  dimensionless  area 
within  C’,  we  finally  obtain 


'B  n  =  l 


(3.16) 


Atlargedistancesdownthechannel.We)appmachesalimitmgminimumvalue. 

If  A.  is  the  smallest  of  the  A’  we  have,  as  r  ^  ,  d,at  4«p  Nu(z)  -  A.  C,2  exp  (-A,  c) 

and  0B  (z)  -  Ci2  exp  (-Aj  z)  so  that 


Nu(oo)  =  ^ 


(3.17) 


3.2.3  Overall  Convective  Heat 
Transfer  Coefficient 


The  average  or  overall  convective  heat  transfer  coefficient  is  defined  through  the 
following  expression: 


Q  =  K{AhbL){AT)u^=^  h,(  Tq-T^)  (Ahbdg) 


(3.18) 


where  ( AT )is  =  log-mean-teraperature  difference  (LMTD). 

The  LMTD  may  be  written  in  terms  of  the  inlet  temperature  Tq  and  the  exit  bulk 


temperature  Ti ,  as  follows: 


{AT  =  (  ^0  T’h,  )  {  Ti  Tw)  ^  1  )  (  Tq—  T^,  ) 


^  [  (  T'o-  r,  )  ] 


Let  us  define  the  dimensionless  LMTD,  as  follows: 


In  ( 


(3.19) 


^LN  = 


Then, 


- 


Tq-  Tk 


(  Obx,  ) 


(3.20) 


(3.21) 


The  expression  for  the  average  convective  heat  transfer  coefficient  can  then  be  written, 
from  Eq.  (3.18): 


^  I  *2  dz 

^LN]  Q 


(3.22) 


In  the  fundamental  equation,  we  have  q  =  hwSi  ATw  Now  q  can  either  be  obtained 
by  integrating  Eq.  (3. 1 3)  from  zero  to  5  or,  alternatively,  it  is  the  heat  given  up  by  the  fluid 
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in  cooling  from  Tq  to  Tm,  so  that  q  =  Awoqc^{To  -  T^).  Equating  these  two  and 
introducing  dimensionless  quantities,  we  have  the  mean  logarithmic  Nusselt  number, 
hLND}j/k,  given  by: 


Nu 


4  d 


LN 


=  -3  <"  (  ) 


(3.23) 


3.3  Summary 

In  this  chapter,  a  series  solution  for  the  mathematical  model  of  temperature 
distribution  in  slip  flow  in  microchannels  has  been  obtained.  Considering  the  given 
boundary  condition,  a  temperature  distribution  in  terms  of  a  generaUzed  Fourier  series 
has  been  derived.  Also,  expressions  for  the  local  and  overaU  Nusselt  numbers  have  been 
obtained.  All  these  expressions  can  be  taken  as  functions  of  the  Knudsen  number  and 
the  Graetz  number.  In  order  to  calculate  either  the  temperature  or  the  Nusselt  numbers, 
the  coefficient  and  A;,  must  be  predetermined. 


CHAPTER  4 


EVALUATION  OF  EIGENVALUES 


In  the  last  chapter,  we  obtained  a  series  solution  for  the  temperature  distribution. 
Also,  expressions  for  the  local  and  overall  Nusselt  numbers  have  been  obtained  as 
functions  of  the  Knudsen  number  and  aspect  ratio.  In  this  chapter,  we  present  a  technique 
for  evaluation  of  eigenvalues  for  the  solution  of  the  heat  transfer  problem  in  sUp-flow, 
since  eigenvalues  must  be  predetermined  for  the  calculation  of  either  the  temperature 
distribution  or  the  Nusselt  numbers.  A  matrix  will  be  constructed  and  a  formulation 
described  to  find  the  coefficients  ^;,.^(m,n).  Based  on  these  i»p,^(m,n)  the  eigenvalues 
can  be  calculated  numerically. 

4.1.  Introduction 

We  consider  the  general  domain  D\  Based  on  the  principle  of  the  method  of 
Galerkin,  let  {(J)ni,n}  be  the  complete  set  of  eigenfunctions  of  the  equation 

With  d(f>  Idv  ’=  -AV>  on  C’.  Any  arbitrary  function  0(x,y)  which  satisfies  these  boundary 
conditions  and  which  possesses  continuous  partial  derivatives  up  to  the  second  order  can 
be  expanded  m  an  absolutely  and  uniformly  convergent  series  in  the  form 

00  00 

.'I'o  (4.2) 
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where 


dx  dy 


^p,q{m,n)  = 


.)=[f 

j  j  D' 


4'm^p,q  dx  dy 


^,q(m,n)  -  0  unless  m=pandn  =  q.Wt  can  make  0  the  solution  of  Eq.  (3.5)  so  that 
multiplying  this  equation  by  and  integrating  over Z)  ’  we  have,  by  Eqs.  (4. 1 )  and  (4.3), 

<5m.«(m,  n)  =  A  j  r(x,  y)0<p^^  dx  dy  (4.5) 

where  r(x,y)  =  vv^y)  /  wp.  If  we  substitute  for  0  by  Eq.  (4.2),  then  Eq.  (4.5)  is  reduced 
to  an  infinite  set  of  homogeneous  linear  algebraic  equations 


w  OD 

p?0  bp^(m,n)} 


ap,q  -  0  (  m,/i  =  0, 1,2,  ...  )  (4.5) 


where 


bpAm 


””^IL 


^  ^Tti^r/P p^q  dx  dy 


m  n.amx  axsooia,=d  with  Eqs.  (4.6)  is  symmauicaJ  since6„(„.„)  = 
elminam  for  a  non-trivial  solulion  gives  an  infinite  detemiinantaJ  equation  /fCA)  =  0 
whose  latent  roots  are  the  eigenvalues  of  Eq.  (3.5).  Dividing  each  row  of  JfA;  by  its 
leading  diagonal  elements,  the  resulung  detemtinant  conveiges  if  the  off-diagonal  sum 
is  absolutely  convergent  and  A  has  no  value  which  makes  a  leading  diagonal  element  zero. 
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If  this  condition  is  satisfied,  then  the  convergence  of  sum  Zp=oSq=o  I  flp,,  I  and, 

moreover,  of  sum  Sp=oXq=o  ^p,q(Pi*i)Ap^  \  Op  q  I  follows.  The  eigenvectors  {apq^"^) 
corresponding  to  a  given  rootil  can  then  be  obtained,  theomUcaUy,  in  terns  of  any 
arbitrary  coefficient;  in  practice,  the  deteimination  of  a  given  eigensolution  is  a  problem 


in  numerical  analysis. 


Eigenvalues  ^  can  detennined  by  the  foUowing  infinite  set  of  homogeneous  linear 
algebraic  equations 

00  00 

^0  Qp^q  =  0  (  m,/j  =  0.1,2,  ...  )  (4.6) 


where 


^,q{m,n)  ~  j  j 


^m^p,q  dx  dy 


^m,n( yirn,n  —  (  / 4a  )  (  nt^  +  a^fl^). 


4i  ^im-p\M-q\  +  9  +  ^  }  (4.10) 


j  oj  ^  ^  ^os(^y)  dx  dy  (ij=0.1,2, ... )  (4.11) 


bp^(m,n)  { 


5^  0  for  all  p,q,m,n  ^  0 

=  0  for  any  =  0 


(4.12) 
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4.2  Expansion  of  Eq.  (4.6) 

Considering  the  condition  Eqs.  (4.8)  and  (4.12)  and  underlining  the  only  one  with 
both  p  =  ra  and  q  =  n,  we  have  the  following  expansion: 

1)  ra  =  0,  n  =  0 

(^iXDAim-Xbo.o)ao,o+(6o,iAo.o-Xbo.i)ao.]+(6o,2Ao.o-Xbo.2)ao.2+(6o.3Ao,o->.bo,3)ao,3+... 
(6i,oAo,o-Xbi^o)ai,o+(6i,iAo,o-^bij)ai,i+(6i,2Ao,o-A.bi,2)ai,2+(6i,3Ao,o-A.bi3)ai,3+... 
(62.oAo.o->vb2,o)a2.o+(62,iAo,o->^b2.i)a2.i+(62,2Ao.o-Xb2,2)a2.2+(62,3Ao.o-5^b2,3)a2,3+... 
(63.oAo,o-Xb3.o)a3.o+(63.iAo,o-Xb3,i)a3.i+(63.2Ao.(^-Xb3.2)a3.2+(63.3Ao.o-Xb3,3)a3,3+... 
=  0 

the  underlined  term  may  not  be  equal  to  zero.  All  terms  of  6p,qAo.o  and  coefficients  of  X, 
bp,q  but  ^,oAn,o  are  zeros  according  to  the  condition  Eqs.  (4.8)  and  (4.12);  therefore, 
(So,oAo,o)ao,o  =  0 


2)  m  =  0,  n  =  1 

(6o.oAo,i-Xbo.o)ao,o+(^£uA£u-Xbo.i)ao.i+(6o.2Ao.i-Xbo.2)ao.2+(6o.3Ao.j-Xbo.3)ao.3+... 
(6i,oAo.]-Xbi.o)ai,o+(6i.iAo.,-Xbi.i)ai,,+(6i.2Ao.i-Xbi,2)ai,2+(6i,3Ao,i-Xbi,3)ai,3+... 
(62,oAo.l-Xb2.o)a2,0+(62.lAo.,-Xb2.i)a2.i+(62.2Ao.,-Xb2.2)a2.2+(62.3Ao.l-Xb2.3)a2.3+... 
(63,oAo,]-Xb3.o)a3.o+(63jAo,i-Xb3.i)a3.]+(63,2Ao.i-Xb3,2)a3,2+(63.3Ao,i-Xb3,3)a3.3+... 
=  0 

therefore. 

(So,iAo,i)aoj  =0 

Similarly, 

3)  m  =  0,  n  =  2 
(6o,2Ao,2)ao,2  =  0 


4)  m  =  0.  n  =  3 


(So,3Ao,3)ao,3  =  0 
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5)  ra  =  1,  n  =  0 

(So.oAi,o-^bo.o)ao,o+(6o.iAi,o-Xbo,i)ao,i+(6o.2Ai,o-Xbo,2)ao.2+(6o.3Ai,o-A.bo,3)ao.3+... 
(^i4iAi^Xbi,o)ai,o+(6i,]Ai,o-^bi,i)ai,i+(6i,2Ai,o-Xbi,2)ai,2+(6i3Ai,o-Xbi,3)ai3+... 
(52,0Al,0-^b2,o)a2.0+(62,iAi.0-Xb2,l)a2,l+(62.2Ai,o-Xb2.2)a2,2+(62,3Ai3-X.b2,3)a2,3+... 
(S3,oAi,o->-b3,o)a3,o+(63,iAi,o-Xb3,i)a3,i+(63,2Ai,o-Xb3,2)a3.2+(633Ai,o-Xb3,3)a33+... 
=  0 

therefore, 

(.5j,nAi,ci  )ai,o=  0 


6)  m  =  1,  n  =  1 

(So.oAi,]-Xbo.o)ao,o+(6o.iAi,i-Xbo.i)ao.i+(6o,2Ai.i-Xbo.2)ao.2+(8o,3Ai.i-Xbo,3)ao,3+... 
(6i.oAi,i-Xbi,o)ai,o+(^i4Ai^-Xbi,i)ai,i+(8i,2Ai.i-Xbi,2)ai.2+(6i,3Ai,i-\bi,3)ai3+... 
(62,oAi,]-Xb2,o)a2,o+(62.iAi,i-Xb2,i)a2.i+(82,2Ai.i-Xb2,2)a2,2+(62,3Ai3-Xb2,3)a2,3+... 
(63,oAi,i-Xb3,o)a3,o+(63jAi,i-Xb3,i)a3,,+(63,2Ai,,-Xb3.2)a3.2+(63.3Ai,i-Xb3,3)a33+... 
=  0 

therefore, 

(^U.Ai,i-Xb  1 , 1  )a  1 , 1 +(-Xb  1 ,2)a  1 ,2+(-Xb  i  ,3)a  i  ,3+.. . 

+(-Xb2,  ]  )a2, 1 +(-Xb2,2)a2,2+(-Xb2,3)a2,3+ ... 

+(-Xb33)a3j+(-Xb33)a33+(-Xb33)a33+ ...  =  0 
or, 

(^i.iAj.i  -Xbl  1  i.])ai.i+(-Xbl  1  i,2)ai.2+(-Xbl li3)ai,3+ ... 

+(-Xbl  1 2,i)a2,]  +(-Xb  1 1 2,2)a2.2+(-Xb  1 1 2.3)a2,3+  ... 

+(-Xbl  l3j)a3j+(-Xbl  l3,2)a3,2+(-Xbl  I3  3)a33+ ...  =  0 


7)  m  =  1,  n  =  2 

(8o,oA,,2-Xbo,o)ao.o+(6o,iAi.2-Xbo.i)ao.i+(6o,2Ai.2-Xbo.2)ao.2+(6o,3Ai.2-Xbo,3)ao.3+... 
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(6i,oAi,2-Xbi,o)ai,o+(6i,]Ai,2->^bi,i)aij+(^x;2Ai^X.bi^2)ai,2+(6i,3Ai,2-^bi,3)ai,3+... 
(S2.oAi.2-Xb2.o)a2,0+(62,lAi.2-?tb2.l)a2,i+(62.2Al.2->^b2.2)a2,2+(62,3Ai,2-Xb2.3)a2.3+... 
(63.oAi,2->.b3,o)a3.o+(63,iAi,2-Xb3,])a3.i+(63.2Ai,2-Xb3,2)a3,2+(83,3Ai,2-Xb3,3)a3,3+... 
=  0 

therefore, 

(-Xbi,i)ai,]+(^i4Ai^-Xbi,2)ai,2+(-Xbi,3)ai,3+... 

+(-Xb2,]  )a2,i+(-Xb2,2)a2,2+(-Xb2,3)a2,3+ ... 
+(-Xb3,i)a3,i+(-Xb3,2)a3,2+(-A.b3,3)a3,3+ ...  =  0 
or, 

(-Xbl2i j)ai,i4-(.$]  7Ai,2  ~^bl2i,2)ai,2+(-Xbl2i,3)ai,3+ ... 

+(-Xb  1 22, 1  )a2, 1 +(-Xb  1 22,2)a2,2+(-Xb  1 22,3)a2,3+ ... 

+(-Xbl23,])a3,]+(->tbl23,2)a3,2+(-Xbl23,3)a3,3+ ...  =  0 


8)  m  =  I,  n  =  3 

(6o.oAi,3-Xbo,o)ao,o+(6o.iAi,3-Xbo.i)ao.i+(6o.2Ai,3-Xbo.2)ao.2+(6o,3Ai.3-Xbo.3)ao.3+... 
(§i,oAi,3-Xbi,o)ai.o+(6i.jAi,3-Xbi.i)ai,,+(6i.2Ai.3-Xbi,2)ai.2+(^i,3Au->^b],3)ai,3+... 
(S2.0Ai,3-Xb2,0)a2.0+(62,lAi,3-Xb2.l)a2.i+(62.2Ai,3-Xb2.2)a2.2+(62.3Aj,3-Xb2,3)a2.3+... 
(63.oAi,3-Xb3.o)a3.o+(63.iAi,3-Xb3.l)a3.i+(63.2Ai,3-Xb3.2)a3,2+(63.3Ai.3-Jlb3.3)a3.3+... 
=  0 

therefore, 

(-Xb  1 ,1  )ai ,  i+(-X.b]  ,2)ai,2+(8j_3Aj^-X.b  i,3)ai  ,3+... 
+(-Xb2,i)a2,]+(-Xb2,2)a2,2+(-Xb2,3)a2,3+ ... 
+(-Xb3,i)a3,]+(-Xb3,2)a3,2+(-Xb3,3)a3,3+ ...  =  0 
or 

(-Xbl3i,i)ai,]+(^,2Ai,2  ~^bl3i,2)ai,2+(-Xbl3i,3)ai,3+ ... 

+(-Xb  1 32, 1  )a2,  ] +(-Xb322,2)a2,2+(-Xb  1 32,3)a2,3+  ... 

+(-Xbl33,])a3,i+(-Xbl33,2)a3,2+(-Xbl33,3)a3,3+ ...  =  0 
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Similarly, 

9)  m  =  2,  n  =  0 

)a2,o=  0 


10)m=:2,n=  1 

(-Xb21]j)ai,i+(-Xb21i,2)a],2+(->.b21i3)ai,3+ ... 
+(^24A24  -^b2l2j)a2,]+(-Xb2l2.2)a2.2+(-Xb2l2,3)a2.3+ ... 
+(-Xb2l3,i)a3j+(-X,b2l3,2)a3,2+(-Xb2l33)a3^+ ...  =  0 


11)  ra  =  2,  n  =  2 

(-X.b22ij)aij+(-Xb22i,2)ai,2+(-Xb22]3)a]^3+ ... 
+(-Xb  1 22.i)a2.]+(^,?A2,2  -Xb  1 22,2)a2,2+(-Xb  1 22,3)32,3+ 
+(-Xb223,i)a3,i+(-Xb223,2)a3,2+(-Xb223,3)a3,3+ ...  =  0 


12)  m  =  2,  n  =  3 

(-Xb23i,i)a],]+(-Xb23i,2)ai,2+(-Xb23i,3)ai,3+ ... 
+(-Xb232,])a2,]+(-Xb232.2)a2.7+(5?,^Ao^  -Xb232,3)a2,3+ ... 
+(-Xb233,,)a3,]+(-Xb233,2)a3,2+(-Xb233,3)a3,3+ ...  =  0 


13)  3  =  2,  n  =  0 

(&,nA3,Q  )a3,o=  0 


14)  ra  =  3,  n  =  1 

(-Xb3 1 1, 1  )ai,i+(-Xb3 1  i,2)ai,2+(-Xb3 1 1,3)31,3+ ... 

+(-Xb3 1 2,1  )a2. 1  +(-Xb3 1 2,2)a2,2+(-Xb3 1 2.3)32,3+ ... 

+(^34A34-^b3l3,i)a3,,+(-Xb3l3,2)a3,2+(-Xb3l3,3)a3,3+  ...  =  0 


15) m  =  3,  n  =  2 

(-X,b32ij)aij+(-)i,b32i^2)^l,2+(-^b32i3)ai3+ ... 

+(-A,b322, 1 )a2, 1  +(-^b322,2)a2,2+(-^b3 22,3)a2.3+  . . . 
+(-X.b323,i)a3.j+(^2;2A3;i-Xb323,2)a3,2+(-Xb323,3)a3,3+ ...  =  0 

16)  m  =  3,  n  =  3 

(-Xb33ij)ai,]+(-Xb33i,2)ai,2+(-Xb33],3)ai,3+ ... 
+(-Xb332,i)a2,i+(-Xb332,2)a2.2+Hb332,3)a2,3+ ... 
+(-Xb333,j)a3,i+(-Xb333,2)a3,2+(S3M3^->^b333,3)a3,3+ ...  =  0 


From  all  these  expansions,  we  can  obtain 

(^Ai,o  )ai,o=  0 


and,  in  matrix  form  and  introducing  =  6nu,An^„  ,  Eq.  (4.6)  becomes 


ei^-Xb\ 1 1 , 

-A/7II12 

-A/7ll,3 

-A/?12,] 

c,,2“A./7l2, 2 

-A/7l2,3 

-AM3i,i 

-A/713,, 2 

^!,,3-A/7l3, 

-^^21,1 

—Xb2 I , 2 

-A^2 1 , 3 

-2/722, , 

-A/722,2 

-A/722,,3 

-a/723,, , 

-a/723,, 2 

-A/72  3,, 3 

-Xbll^i  -^b\222  ->1^1223 

-Xb\^2^  -A1>132.2  ^M32,3 

<?2,I~2Z>212 1  -I^,i?2l2  2  ->1^>212  3 

-Xb222x  ^2,2“^^222_2  ->1^?2223 

-2&232J  -Xb22>22  ^2.7^^232,3 


In  the  next  section,  we  will  deal  with  Eq.  (4.13). 
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4.3  Simplification  of  Eq.  (4.13) 

From  condition  (4.12)  and  bnuip^q  being  non-zero  only  if  m+p  and  n+q  are  both 
even,  Eq.  (4.13)  can  be  further  simplified  as 
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Employing  block  matrices,  Eq.  (4.14)  can  be  further  rewritten  as 


where 


0 

*13 

0 

Q 

*15  •• 

0 

0 

*24 

0  ... 

^31 

0 

^3 

0 

*35  - 

0 

*42 

0 

^4 

0  ... 

*51 

0 

*^53 

0 

its  ••• 

=  0 


(4.15) 


o 

X 

o 

1 

- 1 

X 

0 

X 

0 

■ 

Clrs  1 

0  ^2  0  X  ... 

0  X  0  X  ... 

P^ 

II 

X  0  X3  0  ... 

*/’?  ~ 

X  0  X  0  ... 

Ap  = 

0 

X 

0 

0  X  0  X  ... 

^pA 

m 

. 1 

m 

nJ 

n 

•  •• 

4.4.  Summary 

In  this  chapter,  a  technique  for  calculation  of  the  eigenvalues  occurring  in  the  heat 
transfer  problem  in  slip-flow  has  been  derived  by  constructing  a  matrix.  With  this 
formulation,  any  number  of  eigenvalues  can  be  theoretically  determined.  The  next 
chapter  will  deal  with  the  algorithms  and  the  computational  procedure,  and  the 
calculation  of  the  eigenvalues  will  be  carried  out. 


CHAPTER  5 


COMPUTATIONAL  RESULTS 

In  the  previous  chapter,  the  formulation  for  the  calculation  of  the  coefficients 
occurring  in  the  eigenfunction  was  deteimined.  In  this  chapter,  we  will  discuss  the 
procedures  for  the  evaluation  of  eigenvalues  for  the  heat  transfer  problem  in  slip-flow. 
We  will  calculate  the  eigenvalues  and  discuss  the  computational  results. 

5.1  Procedures  of  Computation 

The  procedures  of  computation  are  as  following; 

1)  Calculate  a\  for  certain  Knudsen  number  Kn  by  Eq.  (2.10)  in  Chapter  2 

a,  mn  a,  =  ^  ,5  j, 

2)  Calculate  coefficients  Ji  j  for  certain  Knudsen  number  Kn  and  aspect  ratio  a  by  Eq. 
(4. 1 1 )  in  Chapter  4 


-  W-.  {-^x)  cos  )  dx  dy 


(ij=0,l,2, ... )  (5.2) 


where 


w/C  =  - 


»  C05  Gi  ip 

22/ - i— 


sin  a-  a  1) 

■  ( - — i - )( 1 _ P _ 

I  +  2  a,.  ^  +  2  Kn  a,-  sin  %■ 
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wo/c  =  -21  -^( — _ X- _ . 

/=ia,  1  +  2 5m2  a.  «  cosh  ^  +  2  Kru  a,  sin 

3)  Deterrnine  bp  q(m,  n)  (or  bmnp  q)  by  Eq.  (4.10)  in  Chapter  4 

bp,gim,n)  =  if  }  (5.5) 

4)  Calculate  Cp^q  =  dp^gAp^q  for  certain  m  and  n  by  Eq.  (4.9) 

\,q(m,n)Ap^q  =  {ji2/4a){m^  +  a^n\  (5 

5)  Determine  truncation  eigenfunction  for  certain  m  and  n  by  Eq.  (4. 1 5)  in  Chapter  4 

0  x„  0  X,.  ...  ir  A,  1 


sink  — 


0 

^13 

0 

XlS  •• 

0 

^2 

0 

X24 

0  .. 

X31 

0 

^3 

0 

X35  •• 

0 

X42 

0 

0  .. 

*51 

0 

^53 

0 

.. 

where 


= 


a:  0  X  0 

0  X  0  X 

X  0  X  0 

0  X  0  X 


Xm/i  — 


'x  0  X  0 

0  X  0  X 
X  0  X  0 

0  X  0  X 


Am  = 


where  x  ore  ( e,.,  -  Xbmn,,,  )  arrd  x  are -Xbnm,^  and  e„  =  for  ceruin  m  and  n. 


45 


For  example,  the  truncated  eigenfunction  for  ra,  n  =  0,1  and  for  m,  n  =  0,1, 2, 3 
respectively 


A](X)=  u:ji=o 


(5.8a) 


A3(X)  —  \Xj\  lAjI—  Ixj/I  IY2I  Lc/jl  =  0  ^5 

and  A  (X)  =  lim  Ap(X)  when  p  ->  00 . 

6)  Evaluate  the  eigenvalue  X  by  the  known  truncated  eigenfunction  such  as  Eq.  (5.8) 

7)  Determine  the  fully-developed  Nusselt  number  by  Eq.  (3.17)  in  Chapter  3 
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5.2  Calculation  of  Cj 

Table  5.1  lists  the  first  eight  values  of  each  ai  for  Kn  from  0.005  to  0.3  by  Eq.(5.1). 

Fig.5.1  shows  thatcivaiy  as  functionsof  Kn.  From  Fig.5.1,  we  can  see  that  the  difference 

of  the  adjacentoi  and  aj+j  becomes  smaUer  as  Kn  increases  butgoes  to  aconstant  Jt  as  both 
Kn  and  i  increase. 


Table  5. 1  The  first  eight  values 


0.005 


Ql 

1.5552 


Q2 

4.6658 


4.6202 


Q3 

7.7764 


7.7012 


Q4 

10.887 


10.783 


as 

13.998 


13.867 


fl6 

17.109 


20.221 


r23.333 


23.128 


22.822 


22.498 


22.348 
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5.3  Calculation  of  d\^ 

IPorKn  —  0.00  and  a  =  1  (for  Kn  =  0.00  and  a  =  1  see  Appendix  B),  and  i,j  =  0,1,2,  ,  10, 

the  (ij  j  are  as  following: 


1.0015  0  -0.27819  0  -0.08309  0  -0.03854  0  -0.02209  0  -0.01431' 
00000000000 
-0.27809  0  0.08285  0  0.02237  0  0.00979  0  0.00542  0  0.00342 

00000000000 
-0.08299  0  0.02237  0  0.00725  0  0.00338  0  0.0192  0  0.00122 

0  0  0  0  0  0  0  0  0  0  0 
-0.03844  0  0.00979  0  0.00338  0  0.01670  0  9.73^-^  0  6.31e-4 

00000000000 
-0.02198  0  0.00542  0  0.00191  0  9.73e-4  0  5.82e-4  0  3.84e-4 

00000000000 
-0.01418  0  0.00343  0  0.00122  0  6.30e-4  0  3.94e-4  0  2.63e-^ 


5.4  Determination  of  bp^g(m,  n)  (or  bmnp^g) 

bp,q{m,  n)  (or  bimip^q)  for  m,n  =  1,2,3  and  p,q  =  0,1,2, ... ,  5  are  as  following: 


0  0  0  0  0 

0.41015  0  -0.06390  0  -0.01428 

0  0  0  0  0 

-0.06390  0  0.01134  0  0.002178 

0  0  0  0  0 

-0.01428  0  0.002178.  0  5.340  r-4 


’0  0  0  0  0  O’ 

0  0  0.3463  0  -0.07818  0 

_  0  0  0  0  0  0 
-  0  0  -0.05256  0  0.01352  0 

0  0  0  0  0  0 

0  0  -0.01211  0  0.002718  0 

0  0  0  0  0  0 

0  -0.06390  0  0.3320  0  -0.08338 

(,,3  0  0  0  0  0  0 

0  0.01134  0  -0.05040  0  0.01425 

0  0  0  0  0  0 

0  0.002178  0  -0.01157.  0  0.002908 


0  0  0  0  O’ 

0  -0.07818  0  0.3268  0 

0  0  0  0  0 
0  0.01352  0  -0.04965  0 

0  0  0  0  0 

0  0.002718  0  -0.01137  0 

0  0  0  0  0  0  ' 

0  -0.01428  0  -0.08338  0  0.3243 

i,l5„=  “  ®  0  0  0  0 

0  0.00218  0  0.01425  0  -0.04933 

0  0  0  0  0  0 

0  5.4  0  0.00291  0  -0.01128 


0  0  0  0  .  0  0 

0  0  0  0  0  0 

0  0.3463  0  -0.05256  0  -0.01211 

0  0  0  0  0  0 

0  -0.07818  0  0.01352  0  0.002718 

0  0  0  0  0  0 


'0  0  0  0  0  O' 

0  0  0  0  0  0 

_  0  0  0.2937  0  -0.06466  0 

0  0  0  0  0  0 
0  0  -0.06466  0  0.01624  0 

0  0  0  0  0  0 


0  0 

0  0 

0  -0.05256 
0  0 

0  0.01352 

0  0 


0  0  0  0 

0  0  0  0 

0  0.2816  0  -0.06914 

0  0  0  0 

0  -0.6194  0  0.01715 

0  0  0  0 


0 

0 

0 

0 

0 

0 


0  0  0  O' 

0  0  0  0 

-0.06466  0  0.2771  0 

0  0  0  0 

0.01624  0  -0.061024  0 

0  0  0  0 


_  0  -0.01207  0  -0.06984  0  0.3126 


0  0.02460 


0.01568  0  -0.06173 


0  0 
0  -0.06373 


0  0.01040 


0.001988 


0  0.3705 


0  -0.05134 


-0.01595. 


0  -0.08289 


0  0.01306 


0.002629 


= 


-0.05333 


0.31922 


-0.06984 


0.012390  0 

0  0 

-0.6294  0 

0  0 

0.01568  0 


0  0.01040 


0  -0.05134 


0  -0.5134 
0  0 
0  0.3076 


0.01306 


-0.06721. 


0  0.01306 


0  -0.06721 


0.01658 
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0  0  0  0' 

0.01239  0  -0.05068  0 

0  0  0  0 

-0.06294  0  0.3033  0 

0  0  0  0 

0.01568  0  -0.06631  0 


blSp^q  = 


0  0 
0  0.001988 
0  0 
0  -0.01160 
0  0 
0  0.002629 


0  0 

0  0.01306 

0  0 

0  -0.06721 

0  0 

0  0.01658 


0  0 

0  -0.05038 
0  0 

0  0.3013 
0  0 

0  -0.06591 


be  equal  to  bnWp^qOT  the  matrices  bmn  should  be  equal  to  bnm'^ .  Using  these 
properties,  we  can  check  the  computational  errors.  For  example,  from  bU^^i  =-0.06390 
and  Z^31i.i  =-0.06373,  we  know  the  third  digit  in  these  values  maybe  not  be  correct 

5.5  Determination  of  Thincation  Eigenfunction 

4.9348-/10.41015  0  A0.06390 

^  0  12.33701-A0.3463  0 

~  **• 

'  A0.06390  0  24.6740 1-A0.3320  ... 


■^0.06373  0  -A0.01040 

0  A0.05333  0 

•>^13  =  •»^3i  =  -A0.01040  0  A0.05134  ... 

If  we  take  the  first  3x3  elements  of  each  matrices  only,  we  have  then 

1^7 1  =  (12.33701-k0.3463)[(4.9348-k0.41015)(24.67401-k0.3320)-X20.063902] 
lA2l  =  (19.73921-X0.3294)[(12.33701-X0.3463)(32.07621-k0.3705)-X20.052562] 
U:jI  =  (32.07621-X0.3705)[(24.67401-X0.3320)(44.41322-X0.3076)-X20.051342] 

kjyl  =  U73I  =  ?t^[(0.06373)(0.05333)(0.05134)  -  (0.0104)2(0.05333)] 

Therefore, 

Aj(X)^\Xi\ 

=  (12.33701-X0.3463)[(4.9348-X0.41015)(24.67401-X0.3320)-X2o.063902]=0 
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and 

A3(X)  =  ]Xi\[X2\\X3\-\x3i\\X2\\Xi3\ 

=  {(12.33701-X0.3463)[(4.9348-X0.41015)(24.67401-X0.3320}-X2o.063902]- 
(19.73921-X0.3294)[(12.33701-X0.3463)(32.07621-X0.3705)-X2o.052562]- 
(32.07621-X0.3705)[(24.67401-X0.3320)(44.41322-X0.3076)-X2o.051342]}- 
{(19.73921-X0.3294)X6[(12.33701-X0.3463)(32.07621-X0.3705)-X2o.052652] 

[(0.06390)(0.05333)(0.05134)  -  (0.0104)2(0.05333)]2}  =  0 

In  the  first  order  approximation  of  eigenfunction  Ai(k),  the  value  of  Xj  can  be 
calculated  approximately  by  the  first  tenn  (4.9348-X0.41015),  because  the  value  of  the 
term  )  can  be  neglected.  It  wiU  be  discussed  in  the  following  section. 

5.6  Evaluation  of  Eigenvalue  A  and 
Fully-Developed  Nusselt  Number 

From  the  above  approximation  of  eigenfunctions,  we  can  evaluate  the  eigenvalues. 
Table  5.2  shows  the  comparison  of  Xj  (a  =  1 )  for  /sTn  =  0.00  with  previously  known  values 
(see  Appendix  C).  From  this  comparison,  we  know  that  the  first  eigenvalues  calculated 
by  the  first  term  are  sufficiently  accurate. 


Table  5.2  Comparison  of  X,(a  =  1)  calculated  by  different  approximation 


Dennis  et  al. 

First  term 

Ai(X) 

A3(X) 

11.91 

12.030528 

11.96289311 

11.962888035 

difference  % 

0 

1.012 

0.44410672 

0.44406406 

Table  5.3  gives  the  computational  results  of  X]  (a)  by  the  first  teim  for  different  aspect 
ratios.  Table  5.4  lists  the  comparison  of  X,  {a)  for  Kn  =  0.00  with  previously  known  values 
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Table  5.3  Eigenvalue  ki(a)  for  different  aspect  ratios 


Kn 

>^1(1) 

}^l(2/3) 

ha/2) 

>t](l/4) 

Xi(l/8) 

0 

12.0305 

12.6418 

13.7316 

17.9174 

22.487 

0.005 

12.219 

12.7896 

13.873 

18.0654 

22.653 

0.01 

12.380 

12.9928 

14.071 

18.2681 

22.839 

0.02 

12.680 

13.1959 

14.266 

18.4708 

23.025 

0.04 

13.204 

13.5795 

14.638 

18.8488 

23.465 

0.06 

13.646 

13.9631 

15.011 

19.2153 

23.905 

0.08 

14.027 

14.3466 

15.383 

19.5876 

24.344 

0.10 

14.359 

14.7302 

15.755 

19.9829 

24.784 

0.20 

15.549 

15.8749 

17.281 

20.8838 

25.467 

0.30 

16.295 

16.6361 

17.618 

21.7846 

26.151 

0.50 

17.189 

17.5510 

18.563 

22.6739 

27.517 

Table  5.4  Comparison  of  X](a)  for  Kn  =  0.00  with  previously  known  values 


h(2/3) 

>^1(1/2) 

Xi(l/4) 

Xi(l/8) 

This  work 

12.0305 

12.6418 

13.7316 

17.9174 

22.487 

Dennis  et  al. 

11.91 

12.49 

13.57 

17.76 

22.38 

difference% 

1.01 

1.21 

1.19 

0.89 

0.48 

and  it  shows  that  all  of  the  differences  are  less  than  1.3  percent  From  this  comparison,  we 
know  that  the  first  eigenvalues  calculated  by  the  first  term  are  sufficiently  accurate  for  all 
different  aspect  ratios. 

Table  5.5  shows  the  computed  results  of  Nuo,(a)  for  different  aspect  ititios  by 
Eq.(5.9),  and  Figure  5.2  shows  Nuoo(a)  as  functions  of  /Cn.  It  shows  that  Nusselt  number 
Nuoo  increase  as  the  Knudsen  number  Kn  increases  and  that  the  values  of  Nu*  increase  as 
the  aspect  ratio,  a,  decreases  for  a  fixed  Kn. 
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Table  5.5  Nuoo(a)  for  different  aspect  ratios 


Kn 

Nu«,(l) 

Nu«(2/3) 

Nuo.(l/2) 

Nu«(l/4) 

Nu«(l/8) 

0 

3.0076 

3.1605 

3.429 

4.4794 

5.6217 

0.005 

3.0548 

3.1974 

3.468 

4.5164 

5.6632 

0.01 

3.0950 

3.2482 

3.518 

4.5670 

5.7098 

0.02 

3.1700 

3.2990 

3.566 

4.6177 

5.7563 

0.04 

3.3010 

3.3949 

3.660 

4.7122 

5.8662 

0.06 

3.4115 

3.4908 

3.753 

4.8038 

5.9761 

0.08 

3.5068 

3.5867 

3.846 

4.8969 

6.0861 

0.10 

3.5898 

3.6826 

3.939 

4.9957 

6.1961 

0.20 

3.8873 

3.9688 

4.320 

5.2209 

6.3668 

0.30 

4.0736 

4.1590 

4.405 

5.4462 

6.5376 

0.50 

4.2970 

4.3878 

4.6408 

5.6685 

6.8792 

Table  5.6  k{a)  for  different  aspect  ratios 


Kn 

^tube 

k(l) 

k(2/3) 

k(l/2) 

~k(l/4) 

k(l/8) 

0 

1 

1 

1 

1 

1 

1 

0.005 

1.015 

1.016 

1.012 

1.011 

1.008 

1.007 

0.01 

1.028 

1.029 

1.028 

1.026 

1.020 

1.016 

0.02 

1.054 

1.054 

1.044 

1.040 

1.031 

1.024 

0.04 

1.099 

1.100 

1.074 

1.067 

1.052 

1.043 

0.06 

1.138 

1.137 

1.105 

1.094 

1.072 

1.063 

0.08 

1.170 

1.169 

1.135 

1.122 

1.093 

1.083 

0.10 

1.198 

1.194 

1.165 

1.149 

1.115 

1.102 

0.20 

- 

1.292 

1.256 

1.260 

1.166 

1.133 

0.30 

— 

1.354 

^  1.316 

1.285 

1.216 

1.163 

0.50 

— 

1.459 

1.388 

1.353 

1.266 

1.224 

Note,  k  —  Nuoo  {Kn)  /Nuoo  (0),  kc  —  for  channel  and  ktuhe  —  for  round  tube 


To  see  the  effect  of  the  Knudsen  number  Kn  on  the  Nusselt  number  Nu«  for  a  certain 
aspect  ratio  with  respect  to  the  case  of  non-slip-flow,  a  ratio  k  =  Nu«o  (Kn)  /Nu«,  (0)  was 
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introduced.  Table  5.6  gives  the  computed  results  of  k{a)  for  different  aspect  ratios. 

Figure  5.3  shows  that  the  ratio  ^  as  a  function  of  the  Nusselt  number  Kn  for  different 
aspect  ratios.  It  shows  that  the  ratio  k  increases  as  Kn  increases  for  all  aspect  ratios  and 
that  the  ratio  k  decreases  as  the  aspect  ratio  decreases  for  0<  Kn  <0.30.  The  maximum 
increase  ofk  (or  the  Nusselt  number  Nuo, )  is  as  follows:  35.4  percent  {Kn  =  0.30)  fora  =  1 ; 

31.6percentfora  =  2/3;28.5percentforo=  l/2;2L6percentforc=  l/4;and  16.3  percent 

for  a  =  1/8.  Figure  5.4  shows  the  comparison  of  Nusselt  number  values  with  that  of  a 
round  tube.  The  result  of  an  aspect  ratio  a  =  1  agrees  with  that  of  the  round  tube  (the 
maximum  difference  is  0.33  percent),  because  the  "aspect  ratio”  of  the  round  tube  can  be 
regarded  as  around  1  while  the  ratio  k  becomes  smaUer  with  a  decreasing  aspect  ratio. 

TTieaspectratiorralso  affects  theNusseltnumberNu„.Fig.5.5  shows  the  ratioof  the 
Nusseltnumberwith  the  aspect  ratio(fl<l)toNusseltnumberwithfl=:l,orNu(a)/Nu(l), 

as  a  function  of  aspect  ratio  a.  As  shown  in  Fig.  5.5,  these  ratios  increase  greaUy  as  the 

aspect  ratio  a  moves  toward  0;  however,  the  effects  are  weaked  slightly  with  the  increase 
of  the  Knudsen  number  Kn. 

5.7  Local  Heat  If-ansfer  Coefficient 
for  the  Case  of  a  Microtube 

The  local  heat  transfer  coefficient  No,  plays  an  important  role  in  detetmining  the 
thermal  effect  of  Kn  in  the  entrance  length.  To  depict  how  the  local  heat  transfer 
coefficient  No,  behaves  with  length,  at  least  two  eigenvalues  ate  needed.  Unfoitunately, 
it  is  extreraly  difficult  at  present  to  obtain  the  second  eigenvalue  in  the  case  of  a 
microchannel.  However,  based  on  the  simUarity  of  the  circular  microtube  and  the  square 
microchannel.  the  situation  in  the  Graetz  problem  with  slip-flow  condition  can  be 
revealed  by  considering  the  microtube  solution. 
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Figure  5.6  shows  the  local  Nux  value  in  the  case  of  a  microtube  as  a  function  of  x*/Gz 
for  Kn  =  0.02  and  with  the  number  of  eigenvalues  as  a  parameter.  The  value  of  the  local 
Nusselt  number  converges  dramatically  with  the  increase  in  the  number  of  eigenvalues  in 
the  computation.  When  x*/Gz  is  0.02,  the  error  in  Nu^  is  0.7  percent  when  two 
eigenvalues  are  used  and,  comparing  to  the  straight  Une  ( using  one  eigenvalue ),  the  error 
is  14  percent.  It  can  be  concluded  that  the  results  using  four  eigenvalues  are  sufficiently 
accurate  for  x*/Gz  >  0.02.  When  x*/Gz  is  greater  than  0.05,  the  error  is  at  most  1.3 
percent  -  that  is,  all  three  plots  become  nearly  flat,  indicating  a  thermally  fully-developed 
condition. 
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Figure  5.7  shows  the  local  Nusselt  number  as  a  function  of  Kn.  It  is  obvious  that  Kn 
has  an  influence  on  the  Nusselt  number.  All  the  plots  in  Fig.  5.7  show  that  the  Nusselt 
number  increases  as  Kn  increases  and  that  this  effect  is  magnified  near  the  entrance.  When 
x*/Gz  is  greater  than  0.05,  all  the  plots  become  nearly  flat,  indicating  a  thermally 
fully-developed  condition. 
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Fig.  5.7  The  local  Nusselt  numbers  as  functions  of  x*/Gz  and  Kn 


From  the  above  discussion,  we  can  see  that: 

( 1 )  slip-flow  has  a  positive  influence  on  the  heat  transfer  coefficient  and  can  enhance 
the  heat  transfer  efficiency; 

(2)  the  influence  depends  on  the  Knudsen  number  and  increases  as  Kn  increases; 

(3)  when  Kn  is  equal  or  greater  than  0.02,  the  increase  in  the  fuUy-developed  Nu  is 
greater  than  5  percent  so  that  the  effect  of  sUp-flow  should  be  taken  into 
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consideration  in  the  computations  of  the  heat  transfer  coefficient;  and 
(4)  that  the  influence  of  Kn  on  Niu  will  decrease  as  Kn  increases. 

5.8  Discussion  on  Other  Eigenvalues 

In  the  expressions  of  JpfA;,  if  we  take  to  be  an  ixi  matrix,  then  there  should  exist 
ip  eigenvalues.  For  example,  for  p  =  1  and  i  =  1,  there  exists  only  one  eigenvalue;  for  p  =  1 
and  i  =  3,  thatis  A  i(A)  in  the  calculations,  there  exists  three  eigenvalues;  and  for  p  =  3  and  i 
=  3,  that  is  A3(X)  in  the  calculations,  there  should  exist  nine  eigenvalues  rather  than  the 
five  shown  in  Table  5.7  (four  eigenvalues  are  missing).  The  calculated  results  are  shown 
in  Table  5.7  (see  Appendix  A  in  details).  From  the  results  in  Table  5.7,  it  is  obvious  that 
the  first  eigenvalue  is  reUable,  while  the  others  are  not  reliable  due  to  truncation  error. 
Therefore,  in  this  research,  only  first  eigenvalues  were  evaluated. 


Table  5.7  Comparison  of  X,(a  =1,  Kn  =  0.00)  calculated  by  different  approximation 


Dennis  et  al. 

AjU)  otXi 

X2 

^3a) 

11.91 

11.962893 

11.962888 

X2 

— 

35.62984 

34.96576 

- 

35.62984 

^3 

71.07 

77.06352 

67.20875 

71.40406 

71.40407 

X4 

— 

- 

121.27253 

113.90901 

113.90902 

^5 

157.9 

— 

- 

180.87005 

121.27253 

5.9  Summary 

In  this  chapter,  the  procedures  for  the  evaluation  of  eigenvalues  for  the  heat  transfer 
problem  through  microchanneU  in  slip-flow  were  developed.  The  first  eigenvalues  and 
the  fully-developed  Nusselt  number  Nu_  were  found  for  different  aspect  ratios.  From  the 
comparison  and  discussion,  it  is  evident  that  the  new  technique  for  evaluation  of  the  first 
eigenvalues  and  the  fully-developed  Nusselt  numbem  of  the  Graetz  Problem  through 
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rectangular  raicrochannels  in  slip-flow  is  coraputionally  effective,  and  the  calculated 
Nusselt  numbers  are  valuable  to  predict  the  heat  transfer  coefficients  for  the 
dimensionless  length  x*/Gz  greater  than  0.05. 


CHAPTER  6 


EXPERIMENTAL  APPARATUS  AND  PROCEDURE 

The  purpose  of  this  work  was  to  detemiine  the  heat  transfer  coefficient  and  Nusselt 
number  in  the  rectangular  raicrochannel  and  microtube  in  a  laminar  flow  region  and 
compare  the  analytical  results  with  experimental  measurements.  Section  6.1  details  the 
test  section  with  the  microchannel  and  microtube.  The  flow  loop  and  the  data  acquisition 
system  are  described.  Section  6.2  discusses  the  design  background  for  the  apparatus. 
Section  6.3  describes  the  procedure  used  to  obtain  the  data  in  the  experiment  It  should 
be  noted  that  while  an  attempt  was  made  to  verify  the  analytical  results  experimentally, 
the  data  obtained  from  the  experiments  was  inconclusive.  It  is  hoped  that  the  details  of 

the  experiments  contained  in  the  next  two  chapters  may  aid  others  what  attempt  similar 
experiments  in  slip-flow. 

6.1.  Experimental  Apparatus 
6.1.1  Test  Section 

The  test  secuon  is  the  part  of  the  experimental  apparatus  where  the  microchannel  and 
microtube  are  tested.  The  microchannel/microtube  were  fabricated  from  3.2  mm  (0.125 
in)  thick  6061-T6  aluminum  bar.  Aluminum  was  chosen  for  its  high  thermal 
conductivity.  Conventional  cutting  and  milling  processes  were  used  to  machine  the 
aluminum  pieces  and  channel  blanks,  to  the  final  25.4  mm  x  1 01 .6  mm  ( 1  in  x  4  in)  shape. 
The  fabrication  technique  was  discussed  in  detail  by  Bailey  (1996). 
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The  microchannel  blank  configuration  is  shown  in  Figure  6.1.  The  macrochannel 
with  dimensions  of  1 17  fim  x  24  |im  x  63.5  mm  was  milled  by  Dr.  Craig  Friedrich  of  the 
Institute  for  Micromanufacturing.  The  milling  machine,  referred  to  as  the  Ultra  Precision 
Milling  Center,  is  a  one-of-a-kind  machine  that  was  specially  built  by  Dover  Instrument 
Corporation  for  the  Institute  for  Micromanufacturing  at  Louisiana  Tech  University. 

The  roughness  of  the  channel  was  obtained  by  using  a  WYKO  Roughness/Step 
Tester  (RST)  along  the  channel  at  30  locations.  The  raw  data  are  shown  in  Appendix  D. 
The  average  roughness  was  shown  to  be  2.28  pm. 

The  WYKO  RST  was  also  used  to  determine  the  channel  width  and  depth.  The 
average  top  width  is  120.43  pm,  and  the  average  bottom  width  is  112.76  pm;  the  average 
depth  is  24.04  pm. 

An  aluminum  cover  was  glued  onto  the  surface  of  the  blank  with  epoxy.  A  70  mm 
long  microtube  (shown  in  Fig.  6.2)  was  glued  into  the  channel  having  dimensions  of  350 
pm  X  350  pm  x  63.5  mm.  Two  pieces  of  lexan  also  were  glued  onto  the  glued  aluminum 
parts  to  form  the  whole  test  section  with  a  water  jacket.  Lexan  was  chosen  based  on  its 
combination  of  the  properties  of  insulation,  strength,  and  machinability.  Care  was 
exercised  to  prevent  the  glue  from  entering  the  microchannel  section. 

6.1.2  Flow  Loop 

Figure  6.3.  shows  the  flow  loop  used  to  conduct  heat  transfer  experiments.  It 
consisted  of  (a)  a  helium  gas  cylinder  ( full  pressure  of 2200  psi ),  (b)  three  flowmeters,  (c) 
two  metering  valves,  (d)  a  heater,  (e)  a  test  section,  (f)  two  vacuum  gages,  (g)  a  vacuum 
pump,  and  (h)  a  constant  temperature  water  bath.  The  test  section,  as  shown  in  Fig.  6.3, 
consisted  of  a  microchannel  and  a  water  jacket  with  two  thermocouples.  As  shown  in  Fig. 
6.4,  two  other  thermocouples  located  at  the  ends  of  the  microchannel  fixed  in  the  test 
fixture  (Bailey,  1996)  were  used  to  measure  the  inlet  and  outlet  conditions  of  the 


a)  the  test  section 


grooves  for  glue  escaping 


Fig.  6.1  Side  cut-away  view  of  the  test  section  and 
the  blank  with  microchannel 

( 117  pm  X  24  pm  X  63.5  mm  ) 


loop  for  micro  heat  transfer  experiments 


microchannel  or  microtube.  The  reading  from  the  thermocouple  in  the  water  jacket  was 
used  to  determine  the  overall  heat  transfer  coefficient. 

6.1.3  Data  Acquisition 

The  data  acquisition  system  consisted  of  a  Gateway  2000  personal  computer 
equipped  with  a  LabVIEW  for  Windows  data  acquisition  system.  The  raw  data  consisted 
of  temperature  readings  (°C)  from  the  thermocouples,  voltage  output  from  the  pressure 
transducers,  and  frequency  from  the  flowmeters.  Converted  data  are  displayed  on  the 
monitor  to  allow  the  user  to  determine  when  a  steady  state  has  been  reached.  The  raw  data 
can  be  sent  to  an  Excel  spreadsheet  when  desired.  Once  in  the  spreadsheet,  the  voltage 
output  can  be  converted  into  pressure,  and  the  frequency  output  can  be  converted  into 
volumetric  flow  rate  by  using  the  calibration  curve  supplied  with  the  meters.  In  addition, 
the  pressure  and  temperature  readings  at  the  exit  of  the  flowmeter  can  be  used  to  convert 
the  volumetric  flow  rate  into  a  mass  flow  rate. 

6.2  Design  Background 

There  are  two  restrictions  that  must  be  achieved  in  the  experimental  apparatus: 

1)  The  flow  must  be  in  the  slip-flow  regime,  so 

0.01<Kn<0.30 

2)  The  flow  must  be  laminar,  so 

Re  <  2300 

6.2.1  Working  Fluid,  Sizing  and  Pressure 

For  a  gas,  the  Knudsen  number  is  found  from: 
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Therefore,  there  are  three  factors  which  may  be  chosen  to  achieve  slip  flow: 

1)  Small  geometry  (small  hydraulic  diameter), 

2)  Low  pressure,  and 

3)  Large  value  of  gas  viscosity. 

Table  6.1  shows  the  corresponding  pressure  ranges  for  different  Dh  for  nitrogen  and 
helium  from  Eq.  (6.1).  Knowing  that  a  raicrochannel  of  10  fim  on  aluminum  could  not 
be  mUled  on  the  Ultra  Precision  Milling  Center,  Dr  =  50  pm  is  chosen.  Thus,  for  the  case 
of  the  microtube,  the  diameter  D  =  50  pm,  and  for  the  case  of  microchannel,  the  dimension 
is  1 00  pm  wide  x  25  pm  deep,  Dh  =  40  pm.  Helium  gas  was  chosen  as  the  working  fluid 
because  its  required  pressure  for  a  given  Knudsen  number  is  three  times  higher  than  that 
of  nitrogen.  Note  that  the  atmospheric  pressure  is  patm  =  14.7  psia  (101.3  kPa).  Thus, 

we  conclude  that  a  vacuum  must  be  used  to  achieve  slip  flow.  Table  6.2  lists  the  fluid 
properties  for  two  gases. 


Table  6.1  Corresponding  pressure  for  different  Dh  for  nitrogen  and  heUum 


Dh:  pm 

10 

50 

100  1 

K 

n 

0.01 

0.3 

0.01 

0.30 

0.01 

0.30 

P 

(psia) 

N2 

9.7 

0.32 

1.94 

0.064 

0.97 

0.032 

He 

^29.1 

0.96 

5.82 

0.190 

^2.91 

0.096 

Table  6.2  Properties  of  nitrogen  and  helium  gases 


gc 

kg-m/N-s^ 

V- 

kg/m  s 

M 

g/mol 

Ru 

J/mol-K 

T 

K 

N2 

1 

1.784x10-5 

28.95 

8.3144 

300 

He 

1 

tT987xlO-5 

4.003 

8.3144 

300 
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6.2^  Flow  Rate 

To  assure  the  gas  flow  in  the  slip  flow  condition,  we  let  pi  =  5.82  psia  (300  torr)  and 
P2  =  0. 19  psia  (9.79  torr)  or  less.  Thus,  we  assume  that  flow  in  such  a  condition  (rarefied 
gas)  should  be  treated  as  a  compressible  flow.  Consequently,  the  estimation  of  pressure 
drop,  Reynolds  number,  and  flow  rate  become  complicated  problems. 

The  exit  pressure  and  inlet  pressure  are  related  by  (Shapiro,  1953): 

-P2-=  (p/p  *h 

Pi  (p/p  *),  (6.2) 

where 

-P-  =  _i  /  Y  +  1 

P*  M  y  2[1  +  (y-1)MV2]  (6.3) 

for  helium  gas  y  =  1.667.  Although  the  flow  can  be  treated  as  choked  flow,  that  is,  M2 
=  1,  Ml  must  be  known  for  the  estimation  of  inlet  pressure,  whUe  the  determination  of 
Ml  (or  the  velocity  of  the  fluid  at  the  inlet)  depends  on  the  inlet  pressure  and  outlet 
pressure.  Also,  the  Reynolds  number  is  a  function  of  pressure  and  may  increase  at  most 
ten  times  of  that  at  STP  (standard  temperature  and  pressure,  300  °K  and  101.325  kPa  or 
1  atm)  (White,  1991). 

The  rough  estimation  of  the  Reynolds  number  is  carried  by  using  the  formulas  for 
incompressible  flow.  The  maximum  mass  flow  rate  to  achieve  laminar  flow  may  be  found 
as  follows: 


Re  = 


-  D  G 


=  2300 


G  =  -22-  =  (2300)(0.01987)(10-^)  _ 
(50)(10-^) 


914.0  kg/s-m2 
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III  =  (914.0)(jt/4)(50  X  10-6)2  =  1.795  x  10-6  kg/s 


The  corresponding  volumetric  flow  rate  at  STP  is: 


^  =  _ra  ^  (1.795)(10-6): 

eo  (0.1625) 


11.04  X  10-6  mVs 


V  —  11.04  cm^/s  =  662.6  cm^/min.  at  STP 


This  IS  the  maximum  flow  rate  for  the  slip  flow  condition,  corresponding  to  a  Reynolds 
number  of 2300.  If  the  flowrate  is  reduced  to  10  cm^/min.  at  STP,  which  is  the  minimum 

flow  rate  for  the  flowmeter  in  the  lab  in  the  Institute  for  Micromanufacturing,  the 
corresponding  Reynolds  number  is: 

Re  =  (2300)(  10/662.4)  =  34.7 

Even  ten  times  this  value  (347)  is  much  less  than  2300,  the  value  of  Reynolds  number  for 

laminar  flow.  Thus,  laminar  flow  should  be  assured  with  the  design  parameters.  Taking 

the  experiment  conducted  by  Yu  ( 1 994)  as  a  reference  ( D,  =  52. 1  fim,  nitrogen,  pj  =  42.64 
psia,  flow  rate  =  9.8 1  mymin.  Re  =  255.91),  the  flow  rate  maybe  higher  than  10  ml/rain. 
Therefore,  the  experiments  should  be  carried  out  without  technical  problems. 

6.3  Experimental  Procedure 

The  test  gas  was  supplied  from  an  industrial  helium  gas  cyUnder  and  the  pressure  of 
the  test  gas  was  reduced  down  to  the  proper  pressure  (about  1 00  psia)  by  the  coarse  valve. 
The  pressure  and  the  flow  rate  were  controlled  precisely  by  the  fine  metering  valve. 

Before  entering  the  test  component,  the  gas  was  heated  to  the  desired  temperature 
range  by  a  heater  strip  with  a  length  of  400  mm  wrapped  around  the  stainless  steel  tube. 
The  mput  power  was  controlled  by  a  temperature  controUer.  In  the  test  section,  water  was 
circulated  on  the  outside  of  the  microchannel  by  a  recirculating  pump  in  the  constant 
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temperature  water  bath. 

The  experimental  procedures  in  the  experiments  were  as  follows: 

1 )  Pump  the  test  section  including  the  steel  pipes  connected  to  the  flow  meter  over  40 
hours  to  check  gas  leakage  and  water  leakage. 

2)  Run  the  acquisition  system  half  an  hour. 

3)  Control  the  metering  valve  to  desired  flow  rate. 

4)  Take  the  data  when  the  flow  reaches  steady-state.  It  took  about  20  minutes  for 
stabilization. 

6.4  Data  Reduction 

6.4.1  Data  Reduction  of  the  MicroChannel 

For  the  heat  transfer  test,  the  heat  transfer  coefficient  of  the  inside  channel  is  derived 
by  the  following  procedures.  From  the  energy  balance  of  the  control  volume  as  shown  in 
Figure  6.5,  we  have 


Fig.6.5  Control  volume  of  heat  transfer  in  microchannel 
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Integrating  Eq.  (6.4)  over  the  entire  length  of  the  channel  yields 


L  ^2-T^2 

^  1^  =  f  dT 
Cp  J  L  j  T-T, 
0 


resulting  in 


^2~^w2  _  .  hAw  . 

T  _T  •  p.  )  (6.6) 

The  wall  temperature  was  considered  constant  along  the  microchannel  so  that  the 
heat  transfer  coefficient  h  was  estimated  by  Eq.  (6.6)  and  the  Nusselt  number  by  Eq.  (6.7) 


Nu  =  -^^ln(I^ 
Awk  Ti-T^2 


Fluid  properties  were  evaluated  at  the  mean  bulk  temperature  of  the  fluid.  The 
Reynolds  number  is  evaluated  as  follows: 


Re  =  Gy 


6.4.2  Data  Reduction  of  the  Microtube 

In  the  case  of  the  microtube,  because  the  tube  is  made  of  polymide  (low  thermal 
conductivity),  the  waU  temperature  will  change  along  the  microtube,  decreasing  from 
Twi  to  Tw2  so  that  the  average  temperature  difference  is  over-estimated  by  Eq.  (6.6). 
Traditionally,  a  correction  factor  is  used  to  correct  the  error,  where  F  is  a  function  of 
parameters  F(Ti,  T2,  Twi,  Tw2)  and  the  function  is  readily  available  from  Holman  (1986). 
The  heat  transfer  process  in  the  tube  involves  a  combination  of  convection  and 

conduction.  The  overall  heat  transfer  is  expressed  in  terms  of  an  overall  heat  transfer 
coefficient  U  as: 
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Q  =  U  Aw  AT 

where  U  is  dependent  on  the  definition  of  AT.  If  AT  =  Tflujd  -  Tw ,  the  h  is  given  by 
j  j  ln(^)  Infe 

- i =  -L  _L  Pq  (6  9^ 

UitDj  hjJiDj  ^  2n:kt  ^  27tkc  ^  ' 

Once  U  (=  h)  is  evaluated  from  Eq.  (6.6),  hj  can  be  calculated  from  the  above 
equation.  Fluid  properties  were  evaluated  at  the  mean  bulk  temperature  of  the  fluid. 


CHAPTER  7 


EXPERIMENTAL  RESULTS 

In  this  chapter,  the  experimental  results  involving  the  microchannel  and  the 
microtube  are  discussed.  It  was  noted  that  there  was  a  significant  discrepancy  between 
the  Nusselt  number  measured  experimentally  and  the  Nusselt  number  calculated  from  the 
analytical  model  involving  slip-flow.  The  fact  that  the  gas  temperature  did  not  increase 
from  inlet  to  outlet  as  heat  was  added  to  the  flow  passage  indicated  that  there  was  some 
problem  in  measurement  of  the  gas  temperature  at  the  exit.  Upon  examination  of  the 
thermocouple  location  at  the  flow  passage  exit,  it  was  noted  that  the  thermocouple 
position  was  not  at  the  centerline  of  the  flow  channel;  therefore,  the  indicated  temperature 
was  not  the  gas  exit  temperature,  in  all  likelihood.  It  is  possible  (but  not  probable)  that 
micrc^specific  mechanisms  caused  the  results  of  the  experimental  study  to  deviate 
significantly  from  the  analytical  model. 

7.1.  Heat  Transfer  with  MicroChannel 

In  this  expenment,  inlet  and  outlet  pressures,  flow  rate,  and  inlet  and  outlet  fluid 
temperatures  were  measured.  These  readings  were  used  to  determine  the  heat  transfer 
coefficient  of  the  microchannel  with  a  certain  Knudsen  number.  It  was  found  that  the 
temperature  drop  of  the  flow  media  depended  on  the  flow  rate,  as  shown  in  Table  7.1. 

In  Figure  7.1,  Nusselt  number  was  plotted  as  a  function  of  Reynolds  number.  Re. 

From  this  experiment,  one  can  observe: 

1 )  Constant  waU  temperature.  As  shown  in  Table  7. 1 ,  the  temperature  differences  of 
wall  are  less  than  0.2°C;  it  is  reasonable  to  regard  the  condition  as  an  isothermal  wall. 
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Nusselt  number,  Nu 
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Table  7.1  Experiment  with  microchannel 


Gv 

ml/min 

"Finiet 

°c 

Xautlet 

°c 

Tw2 

°C 

Nu 

Re 

343.3 

29.90 

31.41 

78.41 

78.56 

80.96 

720 

300.9 

29.62 

33.56 

74.89 

74.96 

43.18 

618 

246.7 

29.84 

32.15 

69.83 

69.88 

20.43 

507 

r  198.6 

31.69 

173.3 

6.59 

356 

149.2 

29.31 

72.13 

72.22 

3.59 

125.9 

259 

101.8 

29.18 

68.69 

208 

*47.3 

29.33 

28.79 

70.65 

0.52 

97 

Fig.  7.1  Nusselt  number  as  a  function  of  Reynolds  number 
(  microchannel  of  117  pm  x  24  pm  x  63.5  mm  ) 
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In  Figure  7.1,  Nusselt  number  was  plotted  as  a  function  of  Reynolds  number,  Re. 

From  this  experiment,  one  can  observe: 

1)  Constant  wall  temperature.  As  shown  in  Table  7.1,  the  temperature  differences  of 
the  wall  are  less  than  0.2°C;  it  is  reasonable  to  regard  the  condition  as  an  isothermal  waU. 

2)  MicroChannel  flow.  Fig.  7.1  shows  that  Nusselt  numbers  were  much  larger  than 
those  of  the  empirical  formula  ( Dittus  et  al.  1930)  for  laminar  flow  if  the  flow  rate  is 
greater  than  200  ml/min  (Re  >  400).  It  is  obvious  that  the  flow  condition  cannot  be 
explained  simply  by  slip-flow  or  even  turbulent  flow.  For  the  case  of  laminar  with  Kn  = 
0.00  and  a  =  1/4,  Nu  =  4.48;  when  a  =  1/8,  Nu  =  5.62 ;  and,  even  in  slip-flow  with  Kn  = 
0.30  and  a  =  1/8,  Nu  =  6.54.  Therefore,  there  is  likely  a  different  mechanism  in  a  flow 
through  a  microchannel/microtube  to  explain  such  a  large  discrepancy.  It  wiU  be 
discussed  in  detail  later  in  the  next  section  of  this  chapter. 

3)  Non-heating-uf^phenomeiion.Theexperunemally  measured  gas  temperature  at 

the  ouUet  was  higher  than  the  inlet  gas  temperature  when  the  flow  rate  was  relatively  large 
(343.3  ml/min.,  for  example).  As  the  gas  flow  rate  was  decreased  to  a  value  less  than  1 25 
ml/min.,  the  gas  tempetature  difference  became  quite  small.  In  fact,  a  negative 
temperature  increase  (i.e.,  a  temperature  decrease)  was  indicated  when  the  flowrate  was 
less  than  47.3  ml/min.  This  phenomenon  was  not  anticipated. 

Upon  examination  of  the  test  apparatus,  it  was  noted  that  the  theimocouple  at  the 

flow  channel  exit  was  not  positioned  exactly  along  the  centerline  of  the  flow  channel.  As 

the  gas  exited  the  test  flow  channel,  instead  of  impinging  directly  on  the  thermocouple, 

the  gas  mixed  with  the  gas  already  within  the  plenum  chamber  at  the  flow  passage  exit. 

The  thermocouple  indicated  the  gas  temperature  in  the  exit  plenum  instead  of  indicating 
the  exit  gas  temperature. 
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This  thermocouple  measurement  error  is  the  most  probable  source  of  the  so-called 
“non-heating-up”  phenomenon  observed  in  the  experiments.  For  high  flow  rates,  the  gas 
in  the  plenum  chamber  was  rapidly  replaced  by  gas  exiting  the  flow  passage;  therefore, 
the  thermocouple  indication  was  more  nearly  that  of  the  exit  gas  temperature.  On  the 
other  hand,  at  low  flow  rates,  the  warm  gas  exiting  from  the  flow  passage  mixed  with  the 
much  larger  volume  of  colder  gas  already  in  the  plenum  chamber.  As  a  result,  the  gas 
mixture  temperature  (measured  by  the  thermocouple)  was  significantly  lower  than  the 
gas  exit  temperature. 

It  is  also  possible  that  there  was  some  leakage  along  the  length  of  the  microchannel, 
although  gas  leakage  would  not  result  in  a  decrease  of  the  gas  temperature  while  the  gas 
was  being  heated. 

7,2,  Heat  Transfer  with  Microtube 

In  this  experiment,  the  flow  rate  was  carefuUy  controlled  to  O.I  ml/min  in  an  attempt 
to  keep  the  inlet  pressure  in  the  vacuum  regime.  After  twenty  minutes,  the  data  were 
measured  with  the  inlet  pressure  737  torr  (14.3  psia)  and  outlet  pressure  5  torr  (0.097 

psia);  the  gas  temperature  in  the  pipe  far  away  from  the  test  section  was  25°C.  The  results 
are  shown  in  Table  7.2. 


Table  7.2  Experiment  with  microtube 


Run 

no. 

Gv 

ral/min 

Tinlet 

°c 

Toutlet 

°c 

Twi 

°C 

Tw2 

°C 

Nu 

Re 

1 

0.11 

29.50 

25.11 

88.41 

82.56 

0.27 

2 

0.11 

30.05 

25.58 

91.35 

85.43 

0.27 

3 

0.11 

30.28 

25.87 

89.50 

84.12 

0.27 

4 

0.11 

30.29 

25.92 

87.73 

82.93 

0.015 

0.27 
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From  Table  7.2,  one  can  observe  that: 

1)  Non-constant  wall  temperature.  It  is  obvious  that  constant  wall  temperature 
condition  cannot  be  achieved  by  the  current  designed  apparatus.  The  temperature 
difference  between  the  inlet  wall  temperature  T^i  and  the  outlet  wall  temperature  Tw2 
(approximately)  is  about  6  °C.  It  is  due  to  that  the  microtube  wall  being  made  of  polymite 
(k  =  0.155  W/m-K)  rather  than  aluminum  (k  =  204  W/m-K).  Therefore,  it  cannot  be 
regarded  as  an  isothermal  wall. 

2)  Non-heating-up-phenomenon.  The  experimental  results  shows  that  the  inlet  gas 
temperature  Tiniet  is  higher  than  the  outlet  gas  temperature  TouUet  by  about  4.5°C,  which 
was  not  the  anticipated  situation.  This  phenomenon  happened  also  in  the  experiment  with 
the  microchannel  when  the  flow  rate  was  47.3  ml/min. 

One  explanation  is  that  the  Joule-Thomson  effect  caused  by  the  expansion  of  the  gas 
at  the  outlet  may  result  in  a  cooling  effect  to  such  an  extent  that  the  outlet  gas  temperature 
Toutiet  is  of  a  smaller  value.  The  expansion  of  a  real  gas  from  a  high  to  a  lower  pressure  in 
an  isenthalpic  (constant-enthalpy)  process  will  result  in  a  temperamre  change,  which 
may  be  either  positive  or  negative,  depending  on  the  pressure  and  temperature  and  type  of 
gas.  For  helium  gas  expanding  from  Tj  =  300  K  =  26.8°C  to  an  exit  pressure  of  1  atm,  we 
find  the  following  values  (Table  7.3)  for  the  exit  temperature,  using  hj  =  h2  (Mann,  1962). 
Note  that  the  temperature  increases  during  the  Joule-Thomson  expansion  for  helium  gas. 
We  may  conclude  that  any  decrease  in  temperature  of  the  helium  gas  could  not  be  caused 
by  the  Joule-Thomson  expansion  effect.  Any  Joule-Thomson  effect  for  helium  gas 
around  the  ambient  temperature  would  result  in  a  warming  effect  or  an  increase  in 
temperature.  Unless  the  inlet  pressure  is  greater  than  about  30  atm  (440  psia),  the 
temperature  change  is  small  (less  than  2°C  or  4°F). 
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Table  7.3  Temperature  change  vs.  pressure  for  helium  gas 


PRESSURE,  pi 
atm  psia  kPa 

EXIT  TEMPERATURE,  T2 

K  °F  X 

(T2- 

“C 

Ti) 

op 

4 

58.8 

405 

300.2 

80.6 

27.0 

6 

88.2 

300.3 

80.8 

27.1 

0.3 

8 

117.6 

811 

300.4 

81.0 

27.2 

10 

147.0 

1013 

300.6 

81.3 

27.4 

0.6 

1.1 

20 

294.0 

2027 

301.2 

82.4 

28.0 

1.2 

2.2 

50 

735.0 

5066 

303.2 

86.0 

30.0 

3.2 

5.8 

Other  explanations  are  as  follows: 

Wang  et  al.  ( 1 994)  pointed  out,  based  on  their  experiments: 

Microscale  heat  transfer  and  transport  phenomena  arc  expected  to  be  quite 
different  from  those  in  customary  situation. 

(1)  For  single-phase  liquid  forced  convection  through  microchannel,  a  fully 
developed  heat  transfer  rcgime  is  initiated  at  about  Re  =  1000-1500.  The 

transition  to  turbulent  mode  is  influenced  by  liquid  temperature,  velocity  and 
micro  size.  . . . 

(2)  Transition  and  laminar  heat  transfer  in  microchannels  are  highly  strange 
and  complicated  [italics  added],  comparcd  with  the  conventionally  sized 
situation.  The  range  of  transition  zone,  and  heat  transfer  characteristics  of  both 

transition  and  laminar  flow  are  highly  affected  by  liquid  temperature,  velocity  and 
micro  size. 

Yu  (1994)  pointed  out  in  his  dissertation: 

One  may  ask  the  reason  and  significance  of  the  shifting.  It  is  well  known  now 
that  the  process  of  convective  heat  transfer  depends  on  flow  field.  The  relative 
rates  of  diffusion  of  heat  and  momentum  are  related  by  the  Prandtl  number.  In  the 
micro  tube  experiments,  it  has  been  shown  that  flow  resistance  was  reduced  both 
in  laminar  and  turbulent  flow  As  it  is  well  known,  if  the  velocity  profile  follows  a 
parabolic  curve,  then  the  friction  factor  f  times  the  Reynolds  number  should  be  64. 
But  in  vanous  micro  tube  experiments,  the  number  is  around  50.  As  discussed 
before  in  Section  1,  it  cannot  be  caused  by  a  slip  boundary  condition.  This 
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manifests  that  there  are  some  alterations  in  velocity  profile,  and  this  in  turn  will 
affect  convective  heat  transfer  [italics  added].  Similarly,  flow  resistance 
reduction  is  also  found  in  turbulent  flow.  All  of  these  will  lead  one  to  expect  that 
heat  transfer  in  micro  tubes  will  behave  differently  both  in  laminar  and 
turbulent  flow  [bolds  added]. 

Also  Bailey  et  al.  (1995)  pointed  out: 


Several  effects  and  conditions  can  exist  in  micro-scale  convection  that  are 
normaUy  neglected  when  considering  macro-scale  flow.  . . .  Another  observed 
micro-scale  effect  is  that  of  large  temperature  variations  of  the  transport  fluid 
[italics  added]  (Wang  and  Peng,  1994).  This  can  cause  a  significant  variation  in 
fluid  properties  throughout  a  micro-system,  invalidating  the  often  used 
assumption  of  constant  properties.  ...  As  of  yet,  it  is  not  completely  clear  when 
these  factors  come  into  play  as  fluid  convection  systems  are  reduced  in  size.  At 
present,  there  are  not  enough  experimental  data  to  make  this  determination. 
Additionally,  there  are  likely  more  micro-specific  effects  and  conditions  that 
have  yet  to  be  observed  [bolds  added]. 

The  experimental  results  for  the  heat  transfer  in  the  microtube  and  the  microchannel 
were  different  in  this  study.  In  both  the  microchannel  and  the  microtube,  the  Reynolds 
number  was  significanUy  less  than  2000;  therefore,  the  flow  pattern  was  completely 
within  the  laminar  flow  regime.  Although  it  is  possible  that  unrecognized  micrc^specific 
phenomena  could  be  responsible  for  the  unusual  experimental  results,  it  is  more  probable 
that  the  placement  error  of  the  thermocouple  at  the  flow  channel  exit  caused  the 
unanticipated  experimental  results  at  the  very  low  flow  rates. 

The  fact  that  the  dimension  of  the  flow  passages  and  the  thermocouple  bead  were  the 
same  order  of  magnitude  could  also  contribute  to  the  thermocouple  measurement  error. 
Because  of  these  uncertainties,  the  experimental  data  obtained  in  this  research  arc  likely 
to  be  of  litde  value  in  assessing  the  validity  of  the  analytical  solutions. 
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7.3  Uncertainty  Analysis  for  the 
Experimental  Data 

Along  with  statistical  analysis,  uncertainly  analysis  mainly  concerns  uncertainty  in 
the  final  results  because  of  uncertainties  in  the  primary  measurements.  The  method 
presented  by  Kline  and  McClintock  (1953)  was  used  in  this  research.  If  the  final  result  R 
is  a  function  of  several  independent  variables,  xi,  X2,  . . .,  Xn, 

R  =  R(xi,X2,  ...,Xn), 

then  the  uncertainty  of  R  associated  with  the  uncertainties  in  measuring  the  primary 
variables  is: 


Ur  = 


or 


From  the  manufacturers’  data,  the  measuring  uncertainties  for  each  independent 
variable  are: 


Ud/D  =  0.0025,  Ul/L  =0.001,  Ucv/Gy  =  0.0128 

Up/P  =  0.0025,  Ut  =  0.4  °C,  Uc/cp  =  U^/p  =  Uv/v  =  Uk/k  =0.005 
(1)  uncertainty  in  the  Reynolds  number 


and  Upe/Re  =  1.49r. 
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(2)  uncertainty  in  dimensionless  temperature 

Ue/0  =  1 .45%  for  the  case  of  microtube;  Ue/S  =  2. 13%  for  the  case  of  amicrochannel  with 
Gv  =  343.3  ml/min  and  Ue/0  =  1.92%  with  Gy  =  149.2  ml/min 
(2)  uncertainty  in  the  Nusselt  number 


0  In  0. 


)2jl/2 


Unu/Nu  =  18.3%  for  the  case  of  a  microtube;  Unu/Nu  =  36.0%  for  the  case  of  a 
microchannel  with  Gy  =  343.3  ml/min  and  Unu/Nu  =  14.9%  with  Gy  =  149.2  ml/min. 


CHAPTER  8 


CONCLUSIONS  AND  RECOMMENDATIONS 

8.1  Conclusions 

In  the  previous  chapters,  the  mathematical  models  of  velocity  distribution  and 
temperature  distribution  were  established,  and  the  expression  for  the  series  solution 
shows  the  importance  of  the  eigenvalues.  Since  those  eigenvalues  were  extremely 
difficult  to  evaluate  direcUy  from  the  original  expansion,  a  concise  matrix  was  derived 
based  on  the  properties  of  the  elements  in  the  original  matrix.  A  truncated  eigenfunction 
was  obtained  which  can  be  used  to  evaluate  the  eigenvalues.  The  procedures  were 
developed  and  some  results  were  obtained.  Also,  the  heat  transfer  experiments  were 
conducted  with  a  single  microchannel  and  with  a  microtube.  From  the  discussions  of 
analytical  and  experimental  results,  the  following  conclusions  can  be  drawn: 

1)  The  technique  for  evaluation  of  the  eigenvalues  of  the  heat  transfer  problem  in 
slip-flow  is  computationally  effective  in  the  evaluation  of  the  first  eigenvalues; 

2)  The  Nusselt  number  increases  as  Kn  increases,  or  the  heat  transfer  is  enhanced 
under  slip-flow  conditions  for  a  given  aspect  ratio; 

3)  When  Kn  is  equal  to  or  greater  than  0.02,  the  increase  in  Nu«  is  greater  than  5 
percent  so  that  the  effect  of  slip  flow  conditions  should  be  taken  into  consideration  in  the 
computations  of  the  heat  transfer  coefficient; 

4)  The  experiments  were  of  no  value  in  assessing  the  validity  of  the  analytical 
solutions  in  this  research. 
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8.2  Recommendations  for  Future  Study 

Considering  the  results  such  as  the  non— heating— up— phenomenon  occurring  in  the 
heat  transfer  experiment  with  a  microchannel/microtube  may  be  caused  totaUy  by 
measurement  errors  due  to  the  small  flow  rate  and  the  fact  that  both  the  dimensions  of  the 
microtube  to  be  measured  and  that  of  the  t-c  bead  are  of  the  same  dimension  level,  it  is 
suggested  that  further  experiments  which  may  result  in  acceptable  data  be  conducted. 
The  suggested  changes  to  the  experiment  include  the  following: 

1)  A  microchannel  machined  in  the  aluminum  should  be  used  to  reach  an 
approximate  isothermal  wall  condition.  The  wall  made  of  a  material  with  high  heat 
conductivity  can  easily  maintain  a  small  temperature  difference  between  the  inlet  and  the 
outlet,  better  approximating  the  isothermal  boundary  condition. 

2)  A  much  shorter  microchannel  or  several  parallel  raicrochannels  should  be 
employed  to  increase  the  flow  rate  for  a  reasonable  pressure  drop. 

3)  A  hot  gas  that  is  cooled  by  a  surrounding  cold  water  jacked  should  be  used  rather 
than  the  reverse  situation.  A  high  temperature  is  desired  at  the  inlet  in  order  to  produce 
a  relatively  large  Kn. 

For  the  analytical  solution,  the  evaluation  of  the  eigenvalues  is  very  important  for 
the  solution  of  the  heat  transfer  problem  through  microchannels  in  slip-flow.  Although 
the  technique  is  effective  for  the  first  values,  it  is  extremely  time-consuming  and  the 
computational  error  is  a  problem  for  other  eignevalues.  It  is  suggested  that  new  codes 
be  developed  on  a  supercomputer  to  solve  the  extremely  large  matrices  to  obtain  more 
eigenvalues  more  effectively  so  that  the  theimal  entrance  effect  in  the  rectangular 
microchannels  can  be  established. 


APPENDIX  A 

CALCULATION  OF  EIGENVALUE  X  FOR  Kn  =  0.00 
WITH  ASPECT  RATIO  a  =  1 
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Calculation  of  Eigenvalue  X  for  Kn  =  0.00  with  Aspect  Ratio  a  =  1  10/27/95 

11  =  P  =  (1+a)/2  and  I2  =  P'  =  (1+1/a)/2 

II  -1  12  =1  a  =1 
i  -0,1..  lOj  =0,1..  10 

Ro  =-i.  y  _j _ 

R(x,y)  =4  y*  — ^ 


1- 


(2  k-  lyn 

..  r 


•tanh 


(2  k-  1)K 


jt-  j^r,  (2  k-  ly  I 


sinj  (2  k-  l)-7t-^  l  i  1  - 


coshl 


I  (2  k-  l)  ii  (x-  0.5) 


I 


a  '  ; 


cosh’^Llll; 
,  2 


Calculating  time  is  about  30  minutes  for  i  =  i  =  10 


•11  ri2 


j 


0  JO 


R(x,y)  /.  x\  !.  y\ 


1  1.00147 

0 

-0.27819 

0 

-0.08309 

0 

-0.03854 

0 

-0.02209 

0 

-0.01431 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-0.27809 

0 

0.08285 

0 

0.02237 

0 

0.00979 

0 

0.00542 

0 

0.00343 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-0.08299 

0 

0.02237 

0 

0.00725 

0 

0.00338 

0 

0.00192 

0 

0.00122 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-0.03844 

0 

0.00979 

0 

0.00338 

0 

0.00167 

0 

9.73085-10"^ 

0 

6.30798-10"'' 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-0.02198 

0 

0.00542 

0 

0.00191 

0 

9.72637- 10"^ 

0 

5.81959-10"^ 

0 

3.83988  10"^ 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-0.01418 

0 

0.00343 

0 

0.00122 

0 

6.29649- 10’^ 

0 

3.93537-10"'' 

0 

2.62711-10"^ 
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® ! 

:  0 

0 

0 

0 

0 

0  j 

’ 

0 

0 

0 

0 

0 

0  I 

'  0 

0 

0 

0 

0 

OJ 

.0 

0 

0 

0 

0 

oJ 

m  =  1  n  =  1 

p  -0,1..5  q  -0.1..5 

u  -  1  J 

D  - id;  -d  4-H  a  \ 

P-<1  4  \  ini-p.:n-q'  im-p  ,n-q  m-rp,n-q  m-i-p,  jn-q/ 


:  0 

0 

0 

0 

0 

0  ■ 

0 

0.41015 

0 

0 

0 

0 

-0.063895 

0 

-0.0142825  ; 

0 

0.41015 

0 

;  0 

0 

0 

0 

0 

0 

0 

0 

0 

:  0 

-0.063895 

0 

0.01134 

0 

0.0021775  i 

0 

-0.063895 

0 

!  0 

0 

0 

0 

0 

0  bll  = 

0 

0 

0 

;  0 

-0.0142825 

0 

0.0021775 

0 

5.4*  10"“  ' 

0 

-0.0142825 

0 

!  0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0  1 

0 

0 

0 

^  0 

0 

0 

0 

0 

i 

0  [ 

.0 

0 

0 

0  0  0  * 

-0.063895  0  -0.0142825  ; 
0  0  0  ' 
0.01134  0  0.0021775  j 

0  0  0^ 

0.0021775  0  5.4  10  '*  ^ 

0  0  0 

0  0  0 

0  0  0 


jO  0  0  0  0  0 

,0  0  0.346255  0  -0.0781775  0 
jO  0  0  0  0  0 

jo  0  -0.052555  0  0.0135175  0 

‘>12  =iO  0  0  0  0  0 

|0  0  -0.012105  0  0.0027175  Oi 

i  0  0  0  0  0  n 


jO  -0.063895  0  0.3319725  0  -0.0833825 
1 0  0  0  0  0  0 

|0  0.01134  0  -  0.0503775  0  0.014245 

‘>13  -jo  0  0  0  0  0 

jo  0.0021775  0  -0.011565  0  0.0029083 
iO  0  0  n  n  A 


u  0  0  0  0  0 

0  0  -0.0781775  0  0.3267675  0 
0  0  0  0  0  0 

0  0  0.0135175  0  -0.04965  0 

bl4  =  0  0  0  0  0  0 

0  0  0.0027175  0  -0.0113742  0  i 

0  0  0  0  A  /\ 


^  0  0  0  0  0 

'0  -  0.0142825  0  -  0.0833825  0  0.324325 
0  0  0  0  0  0 

0  0.0021775  0  0.014245  0  -  0.0493275 

bl5  =!0  0  0  0  0  0 

0  5.410'’  0  0.0029083  0  -  0.0112848 


0  0  0  0  0  0 

0  0  0  0  0  0 

0  0.346255  0  -0.052555  0  -0.012105 

0  0  0  0  0  0 

b21  -  0  -0.0781775  0  0.0135175  0  0.0027175 
0  0  0  0  0  0 

0  0  0  0  0  0 

0  0  0  0  0  0 

-0  0-  0  0  0  0 


0  0  0  0  0  o' 

jO  0  0  0  0  0^ 

10  0  0.2937  0  -  0.06466  0; 

®  0  0  0  oi 

b22  -!o  0  -  0.06466  0  0.016235  0 

0  0  0  0  oi 

'  ®  0  0  0  0  O'' 

0  0  0  0  0 

:0  0  0  0  0  o  ; 
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0 
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0 
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jO 

0 

-0.06466 

0 
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0 

|o 

0 

0 

0 

0 

0  i 

!o 

0 

0 

0 

0 

0 : 

b23 

=  !0 

0.0123898 

0 

-0.0629381 

0 
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io 

0 
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0 
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iO 
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o" 
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0 
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0 

0.0135175 

: 

0 

0  0 

0 

0 

0 
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0 

0 

0 

0 

0 

0  1 

;0  0.3319725 

0 

■0.0503775 

0 

-0.011565  1 

1 

0 

0 
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0 

-0.0619425 

1 

0 

b31  =iO  0 

0 

0 

0 

0 

b32  =  0 

0 

0 

0 

0 

0  1 

.0  -0.083385 

0 

0.0142425 

0 

0.0029107  i 

i 

0 

0 

-0.0691425 

0 

0.0171532 
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0  0  0  0  0  0 

0  0.0104014  0  -  0.051343  0  0.0130546 

0  0  0  0  0  0 

0  -0.0513426  0  0.3075718  0  -0.0672085 
0  0  0  0  0  0  I 

0  0.0130546  0  -  0.0672105  0  0.0165824 

0  0  0  0  0  0 

0  0  0  0  0  0 

0  0  0  0  0  0 


0  0  0  0  0  0  1 

0  0.01134  0  -  0.0503775  0  0.014245  i 

0  0  0  0  0  0 

■0  -  0.0503775  0  0.27003  0  -  0.0662292 

b34  =  0  0  0  0  0  0 

0  0.0142425  0  -  0.0662318  0  0.018148 
0  0  0  0  0  0 

0  0  0  0  0  0 

10  0  0  0  0  0 


|0  0  0  0  0  0 

jo  0.0021775  0  0.014245  0  -0.0493275 

|0  0  0  0  0  0 

.0  -0.011565  0  -0.0662292  0  0.2637127 

b35  =  i  0  0  0  0  0  0 

jO  0.0029107  0  0.018148  0  -0.064789 

1 0  0  0  0  0  0 

!o  0  0  0  0  0 

.0  0  0  0  0  0 
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Comparison  of  different  approximation 


Xl(l)  =(4.9348-  1.0.4I015).(12J370I  -  l..346255).(24.6740]  - 1.. 


used  in  evaluation  of  eigenvalue 


331973)(1)- 


-  (12.33701  -  i,346255).i2.0.0638952 
=10 

«>ot(Xl().), I)  =  11.96289 

>.  =40 

root(Xl(}.).l)  =35.62984 

=70 

n)ot(Xl(l),l)  =77.06352 
^  =30,30.1..  130 


4.9348 

0.41015 

12.33701 

0.346255 

19.73921 

0.2937 


=  12.0317 
=  35.62984 
=  67.20875 


i  =10 

root(X2(i),;)  =34.96576 

J.  =70 

root(X2(l),l)  =67.20875 
i  =120 

roof(X2(l),l)  =  121.27253 


12.33701 

0.346255 

19.73921 

0.2937 

32.07621 

0.281595 


=  35.62984 
=  67.20875 
=  113.90902 


>.  =70,70.05..  190 


X3(X)  =(24.674-  i-0.331973)  (32.07621  -  i  0281595)  (44.41322-  10.270) 
-  (32.07621  -  10.281595)1^  0.0691425^ 


60  80  100 


120  140  150 


1  =70 

n)ot(X3(l),X)  =71.40406 
1  =110 

root(X3(X),l)  =  113.90901 
X  =180 

root(X3(X),X)  =  180.87005 


24.67401 

0.331973 

32.07621 

0.281595 

44.41322 

0.27003 


=  74.32535 


=  113.90902 


=  164.47513 


a  =  h/b  aspect  ratio 


a  =  1  i  =  1 

n  =0,1..5 

m  =0,1..5 


in,n 


4  a 


0 

2.4674011 

9.8696044 

22.2066099 

39.4784176 

61.6850275 


2.467401 1 

4.9348022 

12.3370055 

24.674011 

41.9458187 

64.1524286 


9.8696044 

12.3370055 

19.7392088 

32.0762143 

49.348022 

71.5546319 


22.2066099 

24.674011 

32.0762143 

44.4132198 

61.6850275 

83.8916374 


39.4784176 

41.9458187 

49.348022 

61.6850275 

78.9568352 

101.1634451 


61.6850275 

64.1524286 

71.5546319 

83.8916374 

101.1634451 

123.370055 


APPENDIX  B 

CALCULATION  OF  dy  AND  bmnp^  FOR  Kn  =  0.02 
WITH  ASPECT  RATIO  a  =  l 


Calculation  of  ^  and  bmnp.,  for  Kn  =  0.02  with  Aspect  Ratio  a  =  1 


Kn  -0.02  a  =  .01,. 02 ..22 


-2  Kn  a  tan(a)  -  j 


2000 


•«(“) 


1000  h 


11(a) 


4- 


“-looo 


10 


15 


20 


nx)t(f(a),a)  =  1.510451617057772 
a  =4.533017031227471 
root(f(a),a)  =4.533017098853988 
«  =7.560312429685289 
n)ot(f(a),a)  =7.560312976908389 
a  =10.59472926320076 
n)ot(f(a),a)  =  10.59472964716824 
tt  =  13.63777351248604 
root(f(a),a)  =  13.63777658025087 
a  =  16.69010436230128 

=  16.69010802858231 
a  =19.75169253134426 
n)ol(f(a),a)  =  19.75169677799757 
a  =22.82205687896281 
rool(f(a),a)  =22.82205788967408 


b,  =1-510451617057772 
bj  =4.533017098853988 

bj  =7.560312976908389 

b^  =10.59472964716824 

bj  =  13.63777658025087 

bg  =16.69010802858231 

b,  =  19.75169677799757 
bg  =22.82205788967408 
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wAvO  (w)  for  8  a's  for  ratio  a  =1  (  here  u  =  b )  in  Eq.(19)  [  Eqs.  (12)  and  (13)  ] 


vn  «  sinih.  b. 

wO  =-2-  >  — - '  _[ _ 

i  =  l  (bj)*  l+2  Kn  sin  «  /bj 

*  '  '  ■/  coshl-i 


,'b. 
sinhl  — 


coshi  —  I  i-  2-Kii-h.-stn1i!  — 
\  a  /  *  \  a 


wO  =-0.16265 


wO  --2- 


sini  b.;^ 

\  1/ 


i=l  l  +  lKnsin^bj; 


/b.\ 
sinh!  ~ 

_  \a/ _ 

/b- 

coshi — I  2-Kn-bj'siiihi— - 


i  M,2..8 


wO  =-0.16284 


1  2-Kii-sm^bjj^  ^bjj^ 


/bA  b. 

cosh| — I  4-  2-Kii-bj-siiih  — 


-I  wfO.O.S)  =-0.02681 
wf  1,1)  =-0.00259 
w(0.5,0)  =-0.02498 


11  /  J  - 


Therefore,  using  8  u's  is  accurate  enough 


w(  0,0.5)  =-0.0265 
w(  1,1)  =-0.00305 
w(  0.5,0)  =-0.02629 
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m  =  I  n  =  1 

p  =0,1..5  q  =0,1..5 


b 

1 

id:  ,  1 

-  d 

-  d 

p»q 

4 

\  im-p!  ,  1 

n-q 

im-p 

,n-i~q  in— p,n— q  m-t-p. 

in- 

<1  /' 

'  0 

0 

0 

0 

0 

0 

|0  0 

0' 

0 

0  0 

0 

0.38917 

0 

-0.05173 

0 

-0.01325 

lO  0.38917 

0 

-0.05173 

0  -0.01325 

0 

0 

0 

0 

0 

0 

O 

o 

0 

0 

0  0 

0 

-0.05191 

0 

0.0079 

0 

0.00198 

;0  -0.05191 

0 

0.0079 

0  0.00198 

b  = 

0 

0 

0 

0 

0 

0 

bll  =|o  0 

0 

0 

0  0 

0 

-0.01406 

0 

0.00199 

0 

4.625*10 

-4 

io  -0.01406 

1 

0 

0.00199 

0  4.625  10  '’ 

0 

0 

0 

0 

0 

0 

H 

0 

0 

0 

0  0  i 

0 

0 

0 

0 

0 

0 

iO  0 

0 

0 

0  0  1 

1 

0 

0 

0 

0 

0 

0 

J 

0 

0 

0 

0  0  J 

bj  j  =0.38917 

similarly, 

0 

0 

0 

0 

0 

0’ 

0  0 

0 

0 

0  0 

0 

0 

0.33744 

0 

-0.06498 

0 

0  -0.05173 

0 

0.32419 

0  -0.07135 

0 

0 

0 

0 

0 

0 

0  0 

0 

0 

0  0 

0 

0 

-0.04401 

0 

0.00988 

0 

0  0.0079 

0 

-0.04203 

0  0.01056 

bl2  = 

0 

0 

0 

0 

0 

0 

bl3  = 

0  0 

0 

0 

0  0 

0 

0 

-0.01207 

0 

0.00245 

0 

0  0.00199 

0 

-0.01161 

0  0.00264 

0 

0 

0 

0 

0 

0 

0  0 

0 

0 

0  0 

0 

0 

0 

0 

0 

0 

0  0 

0 

0 

0  0 

.0 

0 

0 

0 

0 

0_ 

.0  0 

0 

0 

0  0 

0 

0 

0 

0 

0 

o' 

iC 

0 

0 

0 

0  0' 

0 

0 

-0.06498 

0 

0.31782 

0 

io  -0.01325 

0 

-0.08323 

0  0.36255 

0 

0 

0 

0 

0 

0 

jo  0 

0 

0 

0  0 

0 

0 

0.00988 

0 

-0.04134 

0 

j  0  0.00198 

0 

0.01306 

0  -  0.05037 

II 

0 

0 

0 

0 

0 

0 

O 

O 

II 

0 

0 

0  0 

0 

0 

0.00245 

0 

-0.01142 

0 

4.625  10'* 

0 

0.00264 

0  -0.01134 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0  0 

0 

0 

0 

0 

0 

0 

1 

io 

0 

0 

0 

0  0 

0 

0 

0 

0 

0 

0_ 

[0 

0 

0 

0 

0  0 
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0 

0 

0 

0 

0 

!0 

1 

0 

0 

0 

0 

0 

iO 

1 

0.33726 

0 

-0.04383 

0 

-0.01127 

iO 

0 

0 

0 

0 

0 

=  io 

-0.06596 

0 

0.00989 

0 

0.00244 

iO 

0 

0 

0 

0 

0 

to 

0 

0 

0 

0 

0 

jo 

0 

0 

0 

0 

0 

■0 

0 

0 

0 

0 

0 

|0 

0 

0 

0 

0 

0 

|0 

0 

0 

0 

0 

0 

jo 

-0.04383 

0 

0.28216 

0 

-0.06079 

0 

I 

0 

0 

0 

0 

0 

=  io 

t 

0.00989 

0 

-0.05363 

0 

0.0132 

10 

! 

0 

0 

0 

0 

0 

|o 

i 

0 

0 

0 

0 

0 

io 

0 

0 

0 

0 

0 

lO 

0 

0 

0 

0 

0 

To 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

io 

t 

-0.01127 

0 

-0.06079 

0 

0.27422 

0 

0 

0 

0 

0 

0 

■=  0 

1 

0.00244 

0 

0.0132 

0 

-0.05239b: 

jo 

0 

0 

0 

0 

0 

10 

1 

0 

0 

0 

0 

0 

iO 

0 

0 

0 

0 

0 

LO 

0 

0 

0 

0 

0 

1 

1 

0 

0 

0 

0 

0 

o' 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0.29343 

0 

-0.0551 

0 

0 

0 

0 

0 

0 

0 

b22  =  0 

0 

-0.05608 

0 

0.01233 

0 

Io 

0 

0 

0 

0 

0 

10 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-0 

0 

0 

0 

0 

0. 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-0.0551 

0 

0.27647 

0 

0 

0 

0 

0 

0 

0 

b24  = 

0 

0 

0.01233 

0 

-0.05277 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.0 

0 

0 

0 

0 

0 

0  0 

0 

0 

0 

0 

0  -0.05191 

0 

0.0079 

0 

0.00198 

0  0 

0 

0 

0 

0 

0  0.32321 

0 

-0.04184 

0 

-0.01081 

0  0 

0 

0 

0 

0 

0  -0.07107 

0 

0.01055 

0 

0.0026 

0  0 

0 

0 

0 

0 

10  0 

0 

0 

0 

0 

!o  0 

0 

0 

0 

0 
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fo 

0  0 

0 

0 

0  ; 

|0 

0 

0  0 

0 

0 

0 

0  -  0.04401 

0 

0.00988 

0  i 

0 

0.0079 

0  -0.04203 

0 

0.01056 

0 

0  0 

0 

0 

1 

0  i 

0 

0 

0  0 

0 

0 

0 

0  0.28136 

0 

-0.05265 

oi 

0 

-0.04184 

0  0.27055 

0 

-0.05815 

0 

0  0 

0 

0 

oi  - 

0 

0 

0  0 

0 

0 

0 

0  -  0.06052 

0 

0.01316 

oi  ! 

0 

0.01055 

0  -  0.05792 

0 

0.0141 

0 

0  0 

0 

0 

0  I  1 

0 

0 

0  0 

0 

0 

0 

0  0 

0 

0 

0 1  i 

0 

0 

0  0 

0 

0 

0 

0  0 

0 

0 

o! 

0 

0  0 

0 

0 

'0 

0 

0 

0 

0 

0 

,0 

1 

0 

0.00988 

0 

-0.04134 

0 

'  0 

0 

0 

0 

0 

0 

iO 

0 

-0.05265 

0 

0.26505 

0 

b34  =  0 

0 

0 

0 

0 

0 

|0 

0 

0.01316 

0 

-0.05697 

0 

jO 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

:0 

0 

0 

0 

0 

0 

10 

i 

0 

0 

0 

0  0 

!o 

0.00198 

0 

0.01056 

0  -0.04105 

lo 

0 

0 

0 

0  0 

1 

-0.01081 

0 

-0.05815 

0  0.26288 

0 

0 

0 

0 

0  0 

0 

0.0026 

0 

0.0141 

0  -0.05657 

0 

0 

0 

0 

0  0 

0 

0 

0 

0 

0  0 

.0 

0 

0 

0 

0  0 

i  =0,1..  10  j  =0,1..  10 


rii  '•12 

M'(x,y) 
j  wO 

•0  Jo  i 


j  iiitcrruptoij^r 


I  0.99883  0  -  0.24674  0  -  0.08526  0  -  0.04342  0  -  0.0222  0  -  0.01381 

:  000000  0  0  0  0  0 

1-0.24538  0  0.06573  0  0.02028  0  0.00912  0  0.00486  0  0.00306 

I  000000  0  0  0  0  0 

1-0.08324  0  0.02024  0  0.00639  0  0.00315  0  0.00162  0  0.00101 

d  =1000000  0  0  0  0  0 

; -0.03838  0  0.00887  0  0.00297  0  0.00158  0  7.94552  10‘‘*  0  4.91934  10''* 

I  000000  0  0  0  0  0 

i -0.02189  0  0.00493  0  0.00169  0  9.45238  10'‘’  0  4.68553  10'‘'  0  2.89189  10"* 

i  000000  0  0  0  0  0 

i-0.01411  0  0.00313  0  0.00108  0  6.24868  10-*  0  3.15992  10-*  0  1.94737-10'* 


m  =0  n  =0 

p  =0,1..5  q  =0,1..5  •*6,0  ‘0  *>7  0  "0  *>8,0  "0 

b  =-  d 

p.q  4  I  Im-pl  ,  In-qi  **1111- pi  ,n^q  ^ _q -  ‘*m4-p.  In- q|  ) 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

o' 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

boo  - 

0 

0 

0 

0 

0 

0  1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

oj 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0  1 

.  0 

0 

0 

0 

0 

0_ 

0 

0 

0 

0 

0 

0  ' 

APPENDIX  C 

THE  FIRST  THREE  EIGENVALUES  FOR  Kn  =  0.00 
BY  DENNIS  ETAL. 
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THE  FIRST  THREE  EIGENVALUES  FOR  Kn  =  0.00  (DENNIS  ET  AL.) 


a 

1.000 

0.667 

0.500 

0.250 

0.125 

h 

11.91 

12.49 

13.57 

17.76 

22.38 

h 

71.09 

51.58 

41.17 

28.17 

25.61 

h 

157.9 

99.71 

94.93 

47.82 

31.81 

0.804 

0.802 

0.789 

0.756 

0.737 

S2 

0.104 

0.064 

0.071 

0.107 

0.091 

0.014 

0.043 

0.020 

0.028 

0.034 

APPENDIX  D 


DATA  OF  ROUGHNESS  AND  DIMENSION  OF  MICROCHANNEL 


I 


9 


DATA  OF  ROUGHNESS  AND  DIMENSION  OF  MICROCHANNEL 


0 


7esf#|  1 

Width 

Depth  1 

1  Top  1  1 

Bottom  1 

1i  115.56 

112 

1  20.97 

r 

133.11 

133.97 


128.01 


118.52 

117.34' 


.2 


118.421 


119.71 
118.41  ■ 


118.52 


117.15 

116.84' 


116.591 


114.69 


114.12 

114.23 


114.971 


114.37 

113.19' 


114:831 


112.6 

112.44 


113.18 


112.41 

112.02 


112 


110.58 


20.58 

22.32 


21.51 ' 


21.96 

22.29 


114.23 

23.67 

112.6 

22.71 

10 

1  117.93 

110.82 

21.55 

11 

119.17 

113.21 

22.85 

12 

122.6 

113.03 

24.61 

13 

128 

112.61 

22.65 

14 

126.12 

113.63 

25.18 

15 

123.2 

113.04 

25.83 

16 

118.89 

112.35 

25.51 

17 

117.22 

111.24 

F  24.86 

24 

114.35 

110.01 

23.271 


24.62 

25.34' 


26.81 


.84 


25.94 


27.221 


23.92 


26.25 


26.66 


24.21 


Standard  G 


Median 


McxJe 


Standard  C 


Variance 


Kurtosis 


Skewness 


ge 


Minimum 


Maximum 


Count  I 


120.4263 


0.991225 


118.52 


5.4291641 


29.47582 


0.910125 


1.396207 


3612.79 


301 


112.7623 


0.22057 


112.605 


112 


1.20811 


1.459529 


-0.0658 


-0.21856 


4.96 


110.01 


114.97 


3382.87 


30 


24.04067 


0.368448 


24.31 


#NI/A 


2.018072 


4.072613 


-0.96811 


-0.01815 


7.23 


20.58 


27.81 


721.22 


30 


APPENDIX  E 


EXPERIMENTAL  HEAT  TRANSFER  DATA 
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Inlet,  outlet  and  wall  temperatures 
Experiment  date:  June  3, 1 996 
Flow  media:  Helium  gas 
Size  of  microchannel:  117  pm  x  24  pm 
Length  of  microchannel:  63.5  mm 
Operation  temperature:  24®C 


Run 

number 

Gv 

ml/min 

Tinlet 

°c 

"^outlet 

°c 

Twi 

°C 

Tw2 

°C 

1 

343.3 

29.90 

37.47 

78.41 

78.56 

2 

300.9 

29.62 

33.56 

74.89 

74.96 

3 

246.7 

29.84 

32.15 

69.83 

69.88 

4 

198.6 

29.77 

31.69 

74.81 

74.89 

5 

173.3 

29.42 

30.54 

70.65 

70.69 

6 

r  149.2 

29.31 

30.12 

72.13 

72.22 

7 

125.9 

29.21 

29.67 

71.07 

71.12 

8 

101.8 

29.11 

29.18 

68.69 

68.74 

9 

47.3 

29.33 

28.79 

70.65 

70.71 
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Inlet,  outlet  and  wall  temperatures 
Experiment  date:  June  18, 1996 
Row  media:  Helium  gas 
Size  of  microtube:  52.1  ftm 
Length  of  microchannel:  70  mm 
Operation  temperature:  25®C 


APPENDIX  F 


REDUCED  HEAT  TRANSFER  DATA 
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Data  reduction  of  heat  transfer 
Experiment  date:  June  3, 1996 
Flow  media:  Helium  gas 
Size  of  microchannel:  1 17  [tm  x  24  pm 
Length  of  microchannel:  63.5  mm 
Operation  temperature:  24®C 


Table  El 


Run 

number 

Gv 

ml/min 

Tinlet 

°c 

Toutlet 

°c 

Twi 

°C 

Tw2 

Nu 

Re 

1 

343.3 

29.90 

37.47 

78.41 

78.56 

80.96 

720 

2 

300.9 

29.62 

33.56 

74.89 

74.96 

43.18 

618 

3 

246.7 

29.84 

32.15 

69.83 

69.88 

20.43 

507 

4 

198.6 

29.77 

31.69 

74.81 

74.89 

11.78 

408 

5 

173.3 

29.42 

30.54 

70.65 

70.69 

6.59 

356 

6 

[  149.2 

29.31 

30.12 

^2.13 

72.22 

i59 

307 

7 

125.9 

r29.21 

29.67 

r  71.07 

71.12 

1.76 

259 

8 

101.8 

29.11 

29.18 

68.69 

68.74 

^  0.11 

208 

9 

47.3 

29.33 

28.79 

70.65 

70.71 

-0.52 

97 

116 


Experiment  date:  June  18, 1996 
Flow  media:  Helium  gas 
Size  of  microtube:  52.1  pm 
Length  of  microchannel:  70  mm 
Operation  temperature:  25°C 

Table  F.2 


APPENDIX  G 

PROPERTIES  OF  THE  MICROCHANNEL  AND  MICROTUBE 
PROPERTIES  OF  HELIUM  GAS 
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Table  G.l.  Physical  properties  of  the  microchannel 


Channel  material 

Aluminum 

Polymide 

Mechanical  Properties 

Density,  kg/ra^ 

2707 

1042 

(Ibm/ft^) 

(169) 

(65) 

103 

Tensile  Strength  (MPa) 

120.7 

(psi) 

(17500) 

(15,000) 

Thermal  Properties 

Thermal  Conductivity 

W/m-K 

204 

0  155 

(Btu/hr-ft-Op)  (118) 

Specific  Heat 

(0.090) 

KJ/kg-K 

0.896 

1  088 

(Btu/lbm 

-®F)  (0.214) 

(0.260) 

Source  of  Data:  J.P.  Holman,  Heat  Transfer,  McGraw-Hill  Book  Co.,  New  York,  1986 


119 


Table  G.2.  Properties  of  heliurn  gas  at  atrnospheric  pressure 
Values  of  p,  k,  and  Cp  are  not  strongly  pressure-dependent  for  He  and  may  be  used  over 
a  fairly  wide  range  of  pressures. 


T,K 

0  , 
kg/m^ 

Cp 

kJ/kgOC 

kg/ra  s 

V 

ra^/s 

k 

W/mOC 

200 

0.2435 

5.200 

15.66x10“^ 

64.38x10-6 

0.1177 

255 

0.1906 

5.200 

18.17 

95.50 

0.1357 

366 

0.13280 

5.200 

1^3.05 

173.6 

0.1691 

477 

0.08282 

5.200 

27.50 

269.3 

0.197 
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